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SUMMARY 

This  document  describes  a method  of  accurately  computing,  in  airborne  fire  control  computers, 
the  impact  range  of  free  falling  projectiles  such  as  bombs  or  bullets.  The  document  has  two  parts. 
In  Part  I older  computation  methods  and  their  limitations  are  described  briefly  and  compared 
with  a new  computation  method,  or  algorithm.  Part  II  presents  a rigorous  mathematical  develop- 
ment of  the  algorithm  with  examples  of  its  accuracy  for  several  typical  Navy  weapons. 

The  new  algorithm  can  contribute  significantly  to  the  flexibility  and  effectiveness  of  Navy  strike 
aircraft.  This  method  of  calculation  is  applied  through  software,  and  it  can  be  used  in  new 
systems  or  reprogrammed  into  existing  systems  without  increasing  hardware  cost. 
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AV  Velocity  deviation  in  one  calculation  interval  (At) 

AY  Altitude  deviation  in  one  calculation  interval  (At) 

0 Dive  Angle 

• Dot  above  symbol  denotes  derivative  with  respect  to  time;  n dots  denotes  n-th  derivative 

p Air  density  (mass) 

pw  Air  density  (weight) 


X 


D 162- 10026-1 


Algorithm 
Ballistic  Wind 

CBU 

CHP 

C'rossrangc 

Dispenser 

Weapon 


Downrange 


Drogue 


I Ejection 


Impact  Range 
RSS 

Runge-Kutta 

Trajectory 

Word 


DEFINITIONS 

A precisely  defined,  step-by-step  method  of  computing  some  quantity. 

A constant,  nonlayered  wind  which  is  equivalent  to  the  actual  wind  in  the  respect 
that  it  affects  the  impact  point  by  the  same  amount. 

Clustered  Bomb  Unit-a  canister-type  bomb  which  dispenses  smaller  bombs  at  some 
point  in  its  trajectory. 

Circular  Error  Probable-the  radius  of  a circle,  centered  about  the  expected  impact 
point,  into  which  a projectile  will  fall  with  a 50*^  probability. 

The  horizontal  distance  measured  in  a direction  90  degrees  clockwise  from  the 
downrange  direction. 

Synonymous  to  CBU. 


The  horizontal  distance  measured  in  a direction  away  from  the  point  of  projectile 
release  and  in  the  vertical  plane  which  includes  the  release  velocity. 

A parachute,  vanes,  or  other  device  deployed  at  the  rear  of  a bomb  to  increase  its 
drag  and  slow  it  down. 

Frequently  bombs  are  physically  ejected  from  an  aircraft  to  ensure  clean  separation. 
The  velocity  imparted  to  a bullet  can  be  considered  an  ejection  velocity. 

The  horizontal  distance  between  the  release  and  impact  of  a projectile. 

Root  Sum  Square- the  square  root  of  the  sum  of  the  squares;  a method  used  to  add 
probabilistic  effects  of  independent  error  sources. 

The  name  of  a class  of  formulas  used  to  evaluate  integrals.  They  are  widely  used  in 
calculating  solutions  of  differential  equations  for  which  analytical  solutions  cannot 
be  found. 

The  path  of  a projectile. 

The  quantity  of  a computer  memory  normally  required  to  store  a single  number  or 
instruction. 
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1.0  INTRODUCTION 

1.1  The  Weapon  Delivery  Problem 

Most  aircraft  weapons  can  be  described  as  ballistic  projectiles.  That  is.  the  only  forces  acting  on 
them  in  flight  are  gravity  and  aerodynamic  drag.  Bullets,  streamlined  bombs,  drogued  bombs, 
clustered  munitions,  and  unguided  rockets  (after  burnout)  are  all  ballistic  projectiles,  (iuided 
weapons  and  weapons  developing  lift  are  not  ballistic  projectiles. 

For  successful  weapon  delivery,  the  pilot  of  an  airplane  must  maneuver  his  plane  into  a position  from 
which  the  weapon  will  fall  to  the  target.  The  problem  is  made  difficult  in  two  ways.  First,  the 
release  must  be  done  accurately  because  the  impact  range  is  very  sensitive  to  release  errors.  Figure  1 
shows  how  the  distance  a bomb  travels  changes  as  bomb  release  conditions  change.  Heading,  dive 
angle,  speed,  altitude,  and  range  from  the  target  must  all  be  exactly  right  when  the  weapon  is 
released.  The  second  difficulty  is  selecting  a correct  combination  of  all  these  factors.  The  path  of  a 
ballistic  projectile  is  curved;  and  the  precise  curvature  depends  on  the  atmospheric  density,  the  size 
and  shape  of  the  projectile,  and  the  speed  and  dive  angle  at  release. 

Several  pilot  aids  and  delivery  tactics  have  been  developed  to  ease  weapon  delivery  problems.  Bomb- 
sights.  release  timers,  and  computers  are  some  of  the  devices;  dive  bombing  is  a delivery  tactic  which 
minimizes  the  effect  of  errors  at  weapon  release. 


Figure  1.  Bomb  Trajectories 
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The  obvious  solution  to  this  problem  is  to  have  the  fire  control  system  calculate  the  correct  range 
to  release  the  weapons  for  any  speed,  altitude,  and  dive  angle.  Because  this  calculation  is  so  difficult, 
no  operational  fire  control  system  can  accurately  compute  range  for  all  the  weapons  and  delivery 
conditions  used.  The  inability  to  perform  this  calculation  accurately  is  a significant  source  of 
weapon  delivery  error  for  all  current  aircraft.  Equally  important,  delivery  tactics  become  constrained 
to  those  for  which  the  fire  control  system  can  do  a fairly  good  job. 

The  subject  of  this  document  is  a new  computing  method,  or  algorithm,  for  the  calculation  of 
weapon  range.  This  method  can  be  used  for  any  ballistic  projectile,  and  it  produces  accurate 
answers  for  any  release  condition  operationally  possible.  This  algorithm  can  be  used  in  current  or 
future  fire  control  systems  having  a digital  computer. 

1.2  History  of  Calculation  Methods 

The  path  of  a ballistic  projectile  can  be  described  by  a relatively  simple-appearing  set  of  differential 
equations.  If  drag  is  left  out,  numerous  complications  can  be  added  to  the  basic  equations  of 
motion  (e.g.,  the  rotation  of  the  earth),  and  a mathematical  solution  can  still  be  found.  This  is  not 
true  once  drag  has  been  included,  because  an  accurate  description  of  the  drag  is  mathematically 
difficult.  Drag  is  a function  of  the  projectile’s  size,  speed,  and  drag  coefficient,  as  well  as  air  density. 
Air  density  varies  primarily  with  altitude  and  secondarily  with  conditions  such  as  the  weather, 
latitude  and  time  of  year.  The  drag  coefficient  varies  with  Mach  number,  which  in  turn  varies  with 
speed  and  air  temperature,  which  varies  with  altitude  and  so  on,  just  as  atmospheric  density  does. 
These  factors  are  empirically  determined  and  are  nonlinear.  At  any  given  speed  or  altitude  the  drag 
acting  on  a projectile  can  be  calculated,  but  an  accurate  function  describing  the  drag  cannot  be  put 
in  the  differential  equations  to  solve  for  the  time  of  fall  or  point  of  impact. 

A way  around  the  lack  of  a solution  has  been  known  for  a long  time.  All  the  bombing  tables  and 
other  ballistic  tables  are  calculated  by  a process  known  as  “numerical  integration.”  In  this  process, 
the  total  path  is  divided  into  many  pieces,  and  each  piece  is  worked  with  an  approximation  that 
all  forces  remain  constant  foi  short  time  periods.  A mathematical  formula  with  constant  forces 
permits  easy  calculation  of  the  position  and  velocity  of  the  projectile  at  the  end  of  0.1  second,  for 
example.  The  forces  on  the  projectile  arc  calculated  for  the  new  point  and  used  for  the  next  0.1 
second.  Figure  2 illustrates  how  this  process  can  be  made  as  accurate  as  desired  by  taking  smaller 
and  smaller  steps. 

A typical  ballistic  table  trajectory  calculation  may  divide  the  fall  of  a bomb  into  a thousand  segments. 
This  yields  very  accurate  answers,  and  each  answer  takes  only  a few  seconds  on  an  IBM  7094  com- 
puter. However,  this  is  too  long  and  the  7094  computer  is  too  large  for  aircraft  fire  control  systems. 
Therefore,  numerous  other  methods  have  been  used  in  the  aircraft  to  provide  weapon  range  informa- 
tion to  the  pilot. 

The  first  advance  in  fire  control  computation  was  the  manually  adjustable  sight  (“depressed  reticle”). 
Before  flight,  a pilot  selects  the  speed,  altitude,  and  dive  angle  he  wants  to  use  to  attack  a target. 

He  goes  to  a ballistic  table  which  gives  the  angle  by  which  the  weapon  will  deviate  from  straight  line 
travel  when  it  reaches  the  ground,  if  released  at  those  conditions.  The  pilot  sets  his  sight  to  this 
angle.  Now  if  he  can  maneuver  the  aircraft  to  the  selected  speed,  altitude,  and  dive  angle  at  the 
instant  the  target  is  in  the  sight,  he  will  get  a hit.  A deviation  in  any  one  of  them,  however,  will 
result  in  a large  miss.  Accurate  weapon  delivery  can  be  done  with  such  a system,  but  it  requires 
excellent  pilot  control.  This  degree  of  control  can  be  attained  only  by  using  intensively  practiced 
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delivery  maneuvers.  Always  using  the  same  maneuver  in  combat  is  dangerous;  moreover,  the  finesse 
achieved  by  a test  pilot  on  a practice  range  is  not  readily  duplicated  by  a green  pilot  under  combat 
conditions.  As  a result,  aircraft  with  relatively  simple  fire  control  systems  may  demonstrate  delivery 
errors  under  100  feet  on  a test  range,  but  produce  average  errors  over  1 ,000  feet  in  combat. 

The  most  modern  systems  (e.g.,  F-1 1 1,  A-6A)  provide  “universal”  mode  weapon  delivery,  using  a 
digital  computer  in  the  fire  control  system  continuously  compute  weapon  range  for  current  condi- 
tions (speed,  altitude,  dive  angle,  wind).  Now  displays  can  be  driven  to  aid  the  pilot  in  steering 
the  aircraft  to  a release  point.  In  one  example,  a heads-up  display  shows  where  the  impact  point  of 
a weapon  would  be  if  released  right  now,  and  how  that  point  will  move  if  the  speed,  heading,  etc., 
are  maintained.  Since  the  pilot  can  also  see  the  target  through  the  display,  he  has  the  information 
needed  to  maneuver  the  airplane  smoothly  to  a release  point,  and  he  maintains  considerable  latitude 
in  just  how  to  do  it.  This  degree  of  flexibility  and  naturalness  contributes  to  a very  significant 
improvement  in  the  operational  usability  of  the  aircraft  and  the  accuracy  of  delivery. 

The  modern  “universal”  systems  are  not  sufficiently  accurate,  however,  and  their  universality  has 
some  severe  limits.  The  inaccuracy  is  due  to  the  previously  discussed  difficulty  of  calculating  a 
ballistic  trajectory.  Approximations  have  been  made  to  keep  the  calculation  simple  and  fast  enough 
for  the  airborne  computer.  Usually  a simplified  form  is  assumed  for  the  drag  function  which  permits 
the  differential  equations  of  motion  to  be  solved  analytically.  This  solution  is  evaluated  by  the 
airborne  computer.  For  a given  type  weapon,  this  simplified  approximation  can  work  well  over  a 
limited  range  of  release  conditions.  Outside  these  conditions,  the  errors  in  the  approximation 
gradually  introduce  larger  and  larger  errors.  For  a different  weapon  or  delivery  regime,  another 
approximation  and  set  of  equations  are  needed.  As  a result,  the  organization  inside  the  fire  control 
computer  begins  to  look  like  older  systems  with  many  “canned”  delivery  modes.  Since  limited 
computer  storage  is  available  for  different  modes,  the  approximations  are  patched  up  to  work  within 
certain  limits  (say  a 100-foot  error)  over  a fairly  wide  range.  Thus  the  approximation  to  truly  i 
universal  operation  is  achieved  at  some  loss  in  accuracy  even  at  those  conditions  for  which  the 
system  is  optimized. 

In  U)67,  algorithms  submitted  in  proposals  by  leading  fire  control  vendors  were  evaluated  at  I he 
Boeing  Company.  The  best  of  these  was  embodied  in  a computer  program  and  tested  over  a wide 
range  of  conditions,  for  both  streamlined  and  drogued  bombs.  Sample  results  are  shown  in  Tables  1 
and  II.  No  existing  fire  control  system  would  include  exactly  these  equation  errors,  but  they  are 
typical  of  the  1967  state-of-the-art.  Note  that  these  errors  come  from  the  equations  alone,  with  no 
allowance  for  round-off  errors,  sensor  errors  or  steering  errors.  It  is  clear  that  refinement  could 
reduce  the  systematic  errors,  but  not  enough  to  make  them  insignificant. 
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Table  I Typical  Equation  Errors  - Ml  17  Bomb 


FLIGHT  PATH  ANGLE  AT  RELEASE 

45°DIVE 

LEVEL 

45°TOSS 

ALTITUDE 

(FT) 

TRUE 

RANGE  (FT) 

EQUATION 
ERROR  (FT) 

TRUE 

RANGE  (FT) 

EQUATION 
ERROR  (FT' 

TRUE 

RANGE  (FT) 

EQUATION 
ERROR  (FT) 

2,000 

1,885 

-15 

10,864 

-330 

28,559 

-718 

5,000 

4,367 

-37 

16,853 

-399 

31,031 

-905 

10,000 

7,832 

-64 

23,333 

-528 

34,626 

-1,226 

15.000 

10,683 

-71 

28.035 

-621 

37,739 

-1,524 

REFERENCE:  NAVAL  WEAPONS  LABORATORY  BALLISTIC  TABLE  NUMBER  101, 
CONDITIONS:  TARGET  AT  SEA  LEVEL 


NO  WINDS 

NO  EJECTION  VELOCITY 
600  KNOTS  RELEASE  SPEED 
-f  MEANS  CALCULATED  RANGE  IS  TOO  LONG 
- MEANS  CALCULATED  RANGE  IS  TOO  SHORT 


Table  II  Typical  Equation  Errors  - Mk  82S  Bomb 


RELEASE  SPEED 

200  KNOTS 

400  KNOTS 

600  KNOTS 

ALTITUDE 

(FT) 

TRUE 

RANGE  (FT) 

EQUATION 
ERROR  (FT) 

TRUE 

RANGE  (FT) 

EQUATION 
ERROR  (FT) 

TRUE 

RANGE  (FT) 

EQUATION 
ERROR  (FT) 

500 

1,507 

-7 

2,483 

+ 78 

3,111 

+185 

1,000 

1,919 

-11 

3,010 

+99 

3,677 

+329 

2,000 

2,360 

-26 

3,554 

+101 

4,255 

+345- 

5,000 

2,901 

-91 

4,244 

+ 23 

5,007 

+267 

REFERENCE:  NWL  BALLISTIC  TABLE  010 
CONDITIONS:  LEVEL  RELEASE 


TARGET  AT  SEA  LEVEL 
NO  EJECTION  VELOCITY 
NO  WIND 
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2.0  NEW  CALCULATION  METHOD 

2.1  Development 

Research  provided  the  clue  which  led  to  a calculation  ureakthrough.  Work  with  mathematical  series 
expansions  showed  better  results  if  the  equations  were  converted  into  a series  with  range,  instead  of 
time,  as  the  independent  variable.  Applied  to  the  numerical  integration  technique,  this  meant  that 
instead  of  dividing  the  trajectory  into  segments  (e.g.,  of  0.01 -second  duration)  the  trajectory  would 
be  divided  into  segments  (e.g.,  covering  10  feet  of  downrange  travel).  The  first  experiments  showed 
good  accuracy  (errors  of  a few  feet  or  less)  could  be  obtained  with  10  or  fewer  integration  steps, 
giving  about  a 100  to  1 improvement  over  the  number  of  steps  used  to  calculate  the  precise  ballistic 
tables.  Two  years  work  produced  refinements  to  extend  the  range  of  weapons  and  delivery  conditions 
which  could  be  handled  by  the  same  basic  method.  Later,  other  refinements  increased  the  speed  of 
computation.  This  development  process  is  continuing,  but  the  new  algorithm  is  better  in  its  current 
state  than  any  computation  method  now  used  in  an  aircraft  fire  control  computer. 

2.2  Description 

The  new  algorithm  is  accurate  and  flexible  because  most  of  the  calculation  is  performed  in  the  same 
manner  that  ballistic  tables  arc  calculated.  At  first,  exactly  the  same  equations  were  used  to  calculate 
atmospheric  density,  speed  of  sound,  gravity  acceleration,  and  drag  coefficient.  Later  work  has 
developed  new  equations  which  yield  the  same  answers,  but  can  be  computed  faster.  At  no  point  in 
the  basic  equations  (the  models  of  the  physical  world)  has  any  approximation  been  used  that  is  dif- 
ferent from  those  accepted  and  used  in  the  calculation  of  ballistic  tables. 

Figure  3 illustrates  how  the  numerical  integration  process  of  the  algorithm  works  and  the  essential 
differences  between  it  and  the  process  used  in  ballistic  table  calculation.  Figure  3,  Point  A,  shows  a 
defined  release  point  (speed,  altitude,  and  dive  angle).  We  know  everything  necessary  about  the 
release  point.  The  impact  range  is  to  be  calculated. 

The  range  of  the  impact  point  must  be  estimated  first  so  the  size  of  the  calculation  intervals  can  be 
chosen.  The  range  the  bomb  would  go  in  a vacuum  is  the  estimate  used.  This  estimate  only  depends 
on  the  release  conditions  an<l  is  calculated  easily  with  a single  equation.  (See  Figure  3,  Point  B.) 

Figure  3,  Point  C,  shows  the  range  estimate  divided  into  five  equal  integration  steps  (an  arbitrary 
number  of  steps  based  on  experimental  results).  This  division  limits  the  amount  of  calculation 
required.  If  total  calculation  time  is  to  remain  small  and  fixed  (a  requirement  for  airborne  use),  the 
number  of  integration  steps  must  remain  small  and  fixed. 

Integration  is  shown  in  Figure  3,  Point  D.  Starting  at  the  release  point,  the  position,  velocity,  etc., 
of  the  projectile  is  calculated  as  it  reaches  the  range  corresponding  to  the  end  of  the  first  step. 

(See  “1 ,”  Point  D.)  No  trajectory  calculations  are  done  at  any  point  between  the  ends  of  the  steps. 

To  calculate  “ 1,”  begin  with  the  release  point  and  calculate  the  average  value  of  the  forces  acting 
on  the  projectile  between  release  and  “1.”  The  simplest  estimate  would  be  that  the  forces  are  the 
same  throughout  the  interval  as  at  release.  This  estimate  is  good  only  if  the  interval  is  very  short 
or  if  the  forces  change  slowly.  A better  estimate  is  needed  to  use  long,  but  few,  steps.  This  is 
obtained  through  the  Runge-Kutta  technique  of  numerical  integration.  The  one  used  in  the  algorithm 
is  the  “standard,  fourth-order,  Runge-Kutta  technique.” 
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The  average  values  of  the  forces  in  the  interval,  expressed  as  functions  of  range,  are  multiplied  by  the 
length  of  the  step  to  find  a close  approximation  of  “1.”  The  same  practice  is  repeated  to  go  from  “1” 
to  “2”  and  so  on  until  “5”  is  reached. 

At  “5”  the  computer  notices  that  the  calculated  position  is  under  ground.  It  is  not  known  where 
between  “4”  and  “5”  the  projectile  struck  the  ground.  Recall  that  the  original  differential  equations 
of  motion  were  functions  of  time.  These  equations  were  rewritten  to  be  functions  of  range.  Now  a 
new  set  is  used,  with  the  same  equations  rewritten  to  be  functions  of  altitude.  In  altitude,  it  is  easy 
I to  decide  the  required  step  length.  The  step  length  is  just  the  difference  between  the  altitude  of  “5” 

and  the  target  altitude.  (See  Figure  3,  Point  E.)  One  more  integration  step,  using  the  same  integra- 
tion technique  with  slightly  different  equations,  carries  the  calculation  to  precisely  the  target  altitude, 
i This  method  of  locating  the  impact  point  is  an  improvement  over  the  way  it  is  done  in  ballistic  table 

calculations.  It  is  obviously  wasteful,  though,  to  calculate  forward  to  “5”  and  back  to  the  impact 
or  return  to  “4”  and  take  a step  down  in  altitude  to  the  impact  point.  A set  of  logical  tests  was 
developed  to  detect  the  last  point  above  ground  and  initiate  the  change  to  the  altitude  coordinate. 

1 (See  Figure  3,  Point  F.)  These  tests  also  detect  cases  (as  in  drogued  weapons)  where  the  trajectory 

has  become  so  steep  that  it  is  more  accurate  to  take  the  remaining  steps  in  altitude  than  in  range. 

The  basic  integration  method  applies  to  any  ballistic  projectile,  but  some  adaptation  is  required  for 
\ each  type  weapon.  The  adaptations  included  in  the  algorithm  are  described  brielly  below: 

' 2.2.1  Streamlined  Bombs 

The  process  described  fits  streamlined  bombs  exactly.  It  is  only  necessary  to  provide  the  bomb  size 
! and  a table  describing  how  the  drag  coefficient  of  the  particular  bomb  varies  with  Mach  number. 

2.2,2  Bullets 

I After  being  fired,  a bullet  is  really  just  a small,  streamlined  bomb  and  is  handled  as  such.  The 

muzzle  velocity  of  a particular  gun-bullet  combination  is  added  to  the  aircraft  velocity  at  release, 
i as  velocity  imparted  to  a bomb  by  an  ejector  would  be  added.  The  weight  and  drag  coefficient 

I table  must  be  provided  as  well  as  muzzle  velocity,  just  as  for  a bomb. 

I 2.2.3  Drogued  Bombs 

These  bombs  are  more  difficult  to  model  accurately.  Their  common  characteristic  is  that  they  are 
released  in  a relatively  low  drag  configuration.  At  some  time  after  release  they  deploy  vanes  or  a 
I parachute  which  greatly  increases  total  drag.  There  is  no  basic  difficulty  in  calculating  bomb 

1 trajectory  before  the  drogue  is  deployed.  Nor  is  there  a basic  difficulty  in  calculating  the  trajectory 

after  the  drogue  is  deployed.  The  problem  is  to  break  the  calculation  correctly  into  two  parts,  one 
using  each  drag  function. 

This  has  been  worked  out  very  satisfactorily  in  the  algorithm  for  the  Mk  82S  bomb.  The  time 
until  the  drogue  opens  is  about  half  a second.  The  procedure  used  is:  (1)  Calculate  the  time  of 
drogue  deployment  using  the  exact  equations;  (2)  Use  a simple  approximation  (Tayloi's  series 
expansion)  to  find  bomb  position  and  velocity  at  the  time  the  drogue  opens;  and  (3)  Use  the 
basic  algorithm  with  the  correct  drag  to  integrate  the  drogued  trajectory  down  to  the  point  of 
impact. 
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Computation  for  other  drogued  bomb  types,  such  as  the  Mk  43,  will  require  adding  appropriate 
logic  for  the  particular  drogue  deployment  mechanism. 

2.2.4  Bluff  Bombs 

These  are  new  bombs,  such  as  the  BLU-58,  which  may  or  may  not  be  drogued.  Bluff  bombs  arc 
treated  satisfactorily  as  streamlined  or  drogued  bombs,  depending  on  whether  they  arc  drogued. 

2.2.5  Cluster  Bombs 

These  weapons  are  released  as  fairly  large,  low-drag. containers.  At  a predetermined  point  in  their 
fall,  they  usi  ome  mechanism  to  dispense  smaller  weapons.  These  smaller  weapons  are  higher  drag 
than  the  container  and  may  also  develop  lift.  The  fire  control  computer  should  predict  the  impact 
location  of  the  center  of  the  pattern  rather  than  individual  positions  of  the  small  bombs.  The  method 
of  applying  the  algorithm  is  to  calculate  the  trajectory  of  the  container  to  the  point  where  it  dis- 
penses the  small  weapons;  switch  to  a drag  function  which  describes  the  motion  of  the  pattern 
center;  and  integrate  down  to  the  ground.  Once  again,  the  important  factor  is  providing  the  logic 

to  switch  the  drag  function  at  the  right  place. 

2.2.6  Rockets 

Work  is  underway  to  extend  the  algorithm  to  include  ungiiided  rockets  and  tracer  bullets.  Rockets 
present  several  new  problems.  They  change  weight  during  fiight;  they  have  thrust  (at  varying  levels) 
as  well  as  drag;  and  they  slew  around  just  after  firing  because  their  launchers  are  not  lined  up  with 
the  aircraft  direction  of  flight.  It  is  felt  that  adequate  provisions  can  be  made  for  these  problems  in 
the  algorithm. 

2.3  Results 

A characteristic  of  numerical  integration  techniques  is  that  they  can  be  made  as  accurate  as  desired 
by  taking  more  integration  steps.  However,  taking  inore  steps  requires  a more  expensive  fire  control 
computer.  A standard  suggested  by  R.  Seeley  of  Naval  Weapons  Center,  China  Lake,  has  been 
adopted:  The  calculation  error  is  insignificant  if  it  is  less  than  10  feet  or  1 milliradian  of  angle 
(0.06  degree)  as  measured  from  the  release  point.  The  arbitrarily  chosen  standard  of  five  integration 
steps  in  range  (plus  one  in  altitude)  provides  answers  which  are  almost  always  within  this  standard, 
even  for  release  conditions  far  outside  the  usual  precision  delivery  envelopes. 

Before  presenting  the  results  obtained  for  some  typical  weapons,  one  characteristic  of  the  tables 
must  be  explained.  The  differences  between  the  calculated  range  and  the  reference  range  have  been 
rounded  to  the  nearest  foot  as  done  in  the  reference  ballistic  tables.  These  two  rounding  processes 
(in  the  ballistic  tables  and  in  the  error  calculation)  introduce  an  uncertainty  of  about  one  foot  into 
the  errors  given.  In  the  case  of  small  errors,  this  gives  an  erratic  appearance  to  the  numbers.  The 
errors,  measured  in  feet,  are  given  as  positive  if  the  calculated  range  is  longer  than  the  reference 
and  as  negative  if  the  calculation  falls  short.  When  the  errors  arc  converted  to  an  angle  for  compari- 
son with  sighting  errors,  the  absolute  value  is  taken.  The  angles  are  rounded  to  the  nearest  tenth 
of  a milliradian  for  the  same  reason  that  the  linear  errors  arc  rounded  to  the  nearest  foot.  Finally, 
the  percentage  error  of  the  calculated  time  of  fall  is  given.  This  can  be  important  in  operation, 
e.g.,  for  wind  correction. 


12 


0162-10026-1 

I 

i 

I 

Table  111  compares  the  new  algorithm  to  that  discussed  earlier,  using  the  M 1 1 7 streamlined  bomb.  j 

Table  IV  compares  the  two  algorithms  tor  the  Mk  82S  drogued  bomb.  Table  V is  for  a bomb  for  ; 

which  the  calculated  ballistic  tables  cover  a wide  range.  The  algorithm  continues  to  give  respectable  | 

results  up  to  the  table  limits.  Mach  2 at  70,000  feet  is  well  outside  the  operational  envelope  of  any 
system  trying  to  precisely  deliver  conventional,  unguided  weapons.  The  generality  of  the  method 
which  achieves  these  results,  however,  provides  the  capability  to  handle  the  new  weapons  or  new 
tactics  not  planned  for  when  a system  is  built,  but  which  become  operational  requirements. 

Table  VI  shows  the  accuracy  achieved  with  a typical  gun  and  round  combination.  Actually,  a 
single  integration  step  gives  fairly  good  answers. 

These  examples  cover  the  range  of  general  weapon  types  to  which  the  algorithm  has  been  applied. 

There  should  be  no  difficulty  in  obtaining  similar  results  for  other  specific  weapons  of  these 
general  types. 


Table  III  Algorithm  Error  Comparison  — M]J7  Bomb 


RELEASE 
ALTITUDE 
(FTt  AND 
ALGORITHM 

RELEASE  ANGLE 

45°  DIVE 

LEVEL 

45°  TOSS 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

(“.) 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

(%) 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

(%) 

2,000  OLD 

-15 

4.0 

0.4 

-330 

5.6 

2.3 

-718 

1.8 

0.9 

2.000  NEW 

+ 1 

0.3 

0.0 

0 

0.0 

0.0 

-27 

0.1 

0.1 

5.000  OLD 

-37 

4.2 

0.2 

-399 

6.6 

1.3 

-905 

4.7 

1.3 

5.000  NEW 

0 

0.0 

0.0 

-1 

0.0 

0.0 

-38 

0.2 

0.1 

10.000  OLD 

-64 

4.0 

0.1 

-528 

8.4 

0.8 

-1,226 

9.8 

1.2 

10,000  NEW 

+3 

0.2 

0.1 

0 

0.0 

0.0 

-63 

0.5 

0.1 

15.000  OLD 

-71 

3.1 

0.3 

-621 

9.4 

0.7 

-1,524 

14.4 

1.2 

15.000  NEW 

+ 2 

0.1 

0.0 

+ 2 

0.0 

0.0 

-85 

0.8 

0.1 

REFERENCE:  NAVAL  HEATONS  LABORATORY  BALLISTIC  TABLE  NUMBER  101.  JULY  1967 


CONDITIONS:  600  KNOT  RELEASE  SPEED 
TARGET  AT  SEA  LEVEL 
NO  WIND 

NO  EJECTION  VELOCITY 
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Table  IV  Algorithm  Error  Comparison  - Mk  82  Snakeye  Bomb 


RELEASE 
ALTITUDE 
(FT)  AND 
ALGORITHM 

RELEASE  SPEED  (KNOTS) 

200 

400 

600 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

Co) 

RANGE 

ERROR 

(FTj 

181 

TIME 

ERROR 

(°«) 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

(%) 

500  OLD 

-7 

1.4 

0.6 

-F78 

5.9 

1.5 

-F185 

3.7 

1.4 

500  NEW 

+1 

0.2 

0.0 

0 

0.0 

0.0 

0 

0.0 

0.0 

1.000  OLD 

-11 

2.3 

0.3 

-F99 

9.5 

0.9 

-F329 

20.9 

0.8 

1,000  NEW 

0 

0.0 

0.1 

-Fl 

0.1 

0.0 

0 

0.0 

0.0 

2,000  OLD 

-26 

5.5 

0.0 

-FlOl 

11.9 

0.7 

-F345 

29.3 

0.7 

2.000  NEW 

0 

0.0 

0.0 

0 

0.0 

0.0 

-1 

0.0 

0.1 

5,000  OLD 

-91 

13.7 

0.3 

-F23 

2.7 

1.0 

-F267 

26.0 

1.5 

5.000  NEW 

-fl 

0.2 

0.0 

-2 

0.2 

0.1 

-1 

0.1 

0.0 

REFERENCE:  NAVAL  WEAPONS  LABORATORY  BALLISTIC  TABLE  NUMBER  010,  OCTOBER  1964 
CONDITIONS:  LEVEL  RELEASE 

TARGET  AT  SEA  LEVEL 
NO  MIND 

NO  EJECTION  VELOCITY 


Table  V Algorithm  Accuracy  - Mk  76  Bomb 


RELEASE  SPEED  (KNOTS) 

RELEASE 

ALTITUDE 

(FT) 

400 

600 

1.200 

RANGE 

ANGULAR 

TIME 

RANGE 

ANGULAR 

TIME 

RANGE 

ANGULAR 

TIME 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

(FT) 

(MR) 

(%) 

(FT) 

(MR) 

(%) 

(FT) 

(MR) 

Co) 

1.000 

0 

0.0 

0.0 

-1 

0.0 

0.0 

-2 

0.0 

0.0 

5,000 

-2 

0.1 

0.1 

-11 

0.2 

0.1 

-40 

0.5 

0.1 

10,000 

-3 

0.1 

0.0 

-22 

0.5 

0.0 

-152  - 

2.2 

0.1 

20,000 

-19 

0.5 

0.0 

-42 

0.8 

0.0 

-91 

1.2 

0.1 

30,000 

-21 

0.4 

0.0 

-48 

0.8 

0.0 

-27 

0.3 

0.0 

40,000 

-22 

0.4 

0.1 

-15 

0.2 

0.2 

-39 

0.4 

0.3 

50,000 

-25 

0.4 

0.2 

-182 

2.2 

3,1 

-73 

0.6 

0.1 

60,000 

-F24 

0.3 

O.l 

-38 

0.4 

0.2 

-36 

0.3 

1.1 

70.000 

-F  7 

0.1 

0.1 

-79 

0.7 

4.6 

-21 

0.1 

1.8 

REFERENCE:  NAVAL  MEAPONS  LABORATORY  BALLISTIC  TABLE  NUMBER  086.  OCTOBER  1966 
CONDITIONS;  LEVEL  RELEASE 
NO  MIND 

TARGET  AT  SEA  LEVEL 
NO  EJECTION  VELOCITY 
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Table  VI  Algorithm  Accuracy  - M61  Gun  with  MSS  Round 


RELEASE 

ALTITUDE 

(FT) 

RELEASE  ANGLE 

10°OIVE 

15°  DIVE 

20°  DIVE 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

m 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

{%) 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

(%) 

500 

0 

0.0 

0.0 

750 

-1 

0.0 

0.0 

1,000 

-3 

0.1 

0.4* 

0 

0.0 

0.7* 

0 

0.0 

1.0* 

1,250 

-1 

0.1 

0.0 

0 

0.0 

0.0 

1,500 

RCF  ua,  .. 

-2 

0.1 

0.3* 

0 

0.0 

0.5* 

1.750 

1 

1 

■1 

0.1 

0.4* 

OF  0.01  SECOND 

ARMAMENT  MEMORANDUM  REPORT  64-5,  FEBRUARY  1964 
400  KNOTS  RELEASE  SPEED 
NO  WIND 

TARGET  AT  SEA  LEVEL 


I 


j 

I 


•TIME  ERROR 
REFERENCE: 
CONDITIONS: 


L. 
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3.0  OPERATIONAL  SIGNIFICANCE 

3.1  Attack  Aircraft  Application 

The  primary  goal  of  the  algorithm  development  is  to  apply  it  to  aircraft  intended  to  deliver  unguided, 
nonnuclear  weapons.  The  best  estimates  show  that  it  will  fit  easily  into  the  new  fire  control  com- 
puters, such  as  that  in  the  A-7E,  in  place  of  the  calculation  methods  now  used.  More  development  is 
required  to  move  the  algorithm  from  the  status  of  a research  tool  to  an  operational  tool.  When  this  is 
achieved,  however,  it  will  be  possible  to  retrofit  the  algorithm  into  existing  aircraft  just  by  changing  the 
computer  program. 

Preparation  of  this  document  is  an  early  step  in  a contract  with  Naval  Weapons  Center,  China  Lake. 
The  primary  objective  of  the  contract  is  to  do  the  mathematical  development  work  necessary  to 
bring  the  algorithm  to  a point  where  it  could  be  programmed  for  the  A-7E  fire  control  computer. 

' 3.2  Contribution  to  System  Performance 

There  is  some  confusion  and  controversy  over  the  accuracy  of  a fire  control  system.  A given  aircraft, 
i on  a test  range  with  the  pilot’s  option  of  delivery  conditions,  might  be  able  to  deliver  streamlined 

( bombs  with  a 100-foot  circular  error  probable  (CEP).  The  same  airplane  may  record  an  average 

■ target  miss  of  1,000  feet  in  North  Vietnam.  The  amount  by  which  the  algorithm  can  improve  the 

i basic  mechanical  ability  of  the  system  can  be  estimated.  What  this  improvement  will  amount  to  in  a 

j combat  situation  is  difficult  to  establish. 

1 

j An  error  in  the  fire  control  equations  will  result  in  a fixed-bias  error.  Over  a reasonable  range  of 

I conditions,  however,  this  bias  error  will  vary  and  may  be  treated-without  gross  injustice-as  a random 

error.  Figure  4 illustrates  the  contribution  to  total  system  error  which  would  be  made  by  a random 
I equation  error  of  50  feet  or  100  feet.  There  are  two  basic  bars:  one  labeled  1968  technology; 

[ another  labeled  1975  technology.  These  represent  calculated  weapon  delivery  errors  which  include 

! all  the  estimated  sensor  errors  and  pilot  errors  (sighting,  steering)  for  first  quality  systems  of  the 

indicated  time  periods.  The  conditions  chosen  for  evaluation  are  shallow-dive  deliveries  of  stream- 
lined bombs  with  release  about  8,000  feet  from  the  target. 

Extensions  above  these  bars  show  how  errors  are  increased  if  the  fire  control  equations  introduce 
random  errors  averaging  50  feet  or  100  feet  in  size.  The  100-foot  figure  is  probably  the  better  esti- 
mate of  current  calculation  error.  Random  errors  add  according  to  the  square  root  of  the  sum  of 
the  squares.  Also,  the  errors  introduced  by  computation  are  not  the  largest  part  of  the  total.  The 
effect  on  system  effectiveness  is  more  apparent  if  one  considers  that  the  chance  of  killing  a point 
target  is  generally  proportional  to  one  over  the  square  of  the  CEP. 

The  change  in  system  accuracy  produced  by  switching  to  an  algorithm  which  produces  very  small 
errors  is  significant.  This  improvement  comes  almost  free  compared  to  the  cost  of  adding  newly 
developed  avionics.  It  is  believed  that  the  change  in  computation  method  will  result  in  a system 
which  also  is  less  susceptible  to  combat  degradation.  The  improved  flexibility  will  permit  the  use 
of  more  varied  tactics,  which  will  contribute  to  improved  aircraft  survival.  A better  chance  of 
survival  should  improve  pilot  performance- if  only  through  practice.  Also,  the  improved  compu- 
tation and  display  techniques  using  it  can  lighten  the  pilot’s  workload  during  the  critical  weapon 
delivery  run. 
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CALCULATED  TEST  RANGE  ACCURACY  WITH  DELIVERY  TACTICS  TYPICAL  OF  S.E.  ASIA 
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Figure  4.  Contribution  oi  Computation  Errors  to  Systom  Errors 
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PART  II 

TECHNICAL  DESCRIPTION 
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1.0  INTRODUCTION 

Part  H presents  the  mathematical  aspects  of  finding  the  impact  range  of  a projectile.  Major  topics 
covered  in  this  part  are : (1)  the  basic  mathematical  model,  (2)  different  coordinate  systems, 

(3)  integration  process,  (4)  adaptations  to  various  weapons,  and  (5)  sample  results.  Coverage  of 
these  topics  leads  to  an  algorithm  for  calculating  weapon  impact  ranges. 

The  basic  technique  used  in  the  algorithm  is  to  transform  the  differential  equations  so  that  down- 
range  is  the  independent  variable.  Then  the  fourth-order,  Runge-Kutta  integration  formulas  are 
used  to  integrate  these  differential  equations  to  obtain  the  impact  range.  Two  variations  of  the 
basic  technique  are  presented  for  special  cases:  (1)  when  the  weapon  trajectory  becomes  very 
steep,  it  is  desirable  to  switch  the  integration  process  to  a set  of  differential  equations  with  altitude 
as  the  independent  variable;  (2)  when  dealing  with  time-fuzed  dispenser  weapons  it  is  convenient 
to  integrate  the  differential  equations  with  time  as  the  independent  variable  over  the  segment  of 
trajectory  between  release  and  fuzing. 

Application  of  the  algorithm  presented  in  this  document  depends  on  the  availability  of  a current, 
advanced,  airborne,  digital  computer.  Also,  it  depends  somewhat  on  the  fact  that  the  weapon  impact 
range  will  be  calculated  repetitively.  This  proeess  provides  a good  estimate  of  the  integration  interval 
size  and,  consequently,  saves  calculation  time. 
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2.0  MATHEMATICAL  MODELS  AND  THEIR  INTEGRATION 

The  differential  equations  describing  the  motion  of  a projectile  are  discussed  in  this  Section.  These 
equations  are  given  first  with  time  as  the  independent  variable.  Then  they  are  transformed  twice; 
once  to  make  the  downrange,  X,  the  independent  variable  and  again  to  make  the  altitude,  Y,  the 
independent  variable.  Also  discussed  are  the  basic  method  used  for  integration  and  the  comparative 
eHlciencies  resulting  from  the  use  of  the  different  coordinate  systems.  In  this  section,  the  discussion 
is  limited  to  streamlined  bombs.  Adaptations  to  other  weapons  are  given  in  Section  3.0. 

2. 1 Basic  Differential  Equations 

The  development  of  an  effective  weapons  release  system  is  inherently  dependent  upon  obtaining 
solutions  of  the  equations  for  the  motion  of  a projectile  within  the  atmosphere.  This  is  generally 
a difficult  mathematical  problem  which  has  not  been  solved  completely.  The  major  difficulty  stems 
from  the  nonlinearities  introduced  by  the  atmospheric  effects  on  a falling  weapon. 

In  choosing  the  mathematical  model,  two  considerations  have  been  kept  in  mind.  The  major  objec- 
tive of  the  mathematical  analysis  is  to  yield  the  weapon  impact  point.  The  main  effect  of  this  is 
that  the  weapon  mass  can  be  assumed  to  be  a point  mass.  Also,  the  choice  of  the  model  is  dictated 
by  the  need  to  evaluate  results  against  some  standard.  Since  the  armed  forces  publish  range  tables 
for  various  weapons,  the  model  is  chosen  to  conform  as  closely  as  possible  to  the  model  used  tor 
these  tables. 

The  equations  of  motion  are  developed  assuming  the  projectile  is  a point  mass  acted  on  only  by  the 
force  of  gravity  and  the  retardation  forces  due  to  air  resistance.  The  trajectory  can  be  restricted  to 
a plane  by  ignoring  crosstrack  effects  such  as  winds.  For  practical  applications,  the  effect  of  winds 
can  be  accounted  for  in  a straightforward  manner. 

The  assumptions  adopted  are  summarized  below: 

a.  The  Earth  is  flat  and  nonrotating. 

b.  The  gravitational  attraction  is  a function  of  altitude. 

c.  The  projectile  is  a point  mass. 

d.  The  projectile  is  not  powered  and  lias  a constant  mass. 

Under  these  assumptions,  the  differential  equations  of  motion  have  the  following  form.  (See 
Reference  1 for  a complete  discussion.) 


^.11^=0 

dt“  dt 


dt-^  dt 


+ g(Y)  = 0 


(2-1) 
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where  X , Y , and  t denote  downrange,  altitude,  and  time,  respectively.  The  coefficient  H , 
which  is  the  drag  function,  is  given  by 


= (2-2) 

where  P is  the  atmospheric  mass  density,  M is  the  bomb  mass,  D is  the  bomb  diameter,  is 
the  weapon  coefficient  of  drag,  and  V is  the  velocity.  The  atmospheric  density  P is  given  as  a 
function  of  altitude  which  is  fitted  to  measured  values  of  atmospheric  density.  , which  is  equal 
to  , is  empirically  derived  and  given  in  tabular  form  as  a function  of  mach  number.  The 
gravitational  acceleration  of  a point  mass  above  the  surface  of  the  Earth  is  given  by 


g(Y)  = go 


(Re  + Y)- 


wherc  g^  is  the  local  value  of  gravity  at  sea  level  and  Re  is  the  local  value  of  the  Earth’s  effective 
radius. 

The  above  differential  equations  are  not  analytically  integrable  if  an  accurate  model  of  H is  used, 
because  an  accurate  model  would  render  them  extremely  nonlinear. 

I he  two  second-order,  differential  equations  given  in  (2-1 ) will  now  be  rewritten  as  four  first-order, 
differential  equations.  This  is  done  to  get  the  differential  equations  in  a form  which  is  more 
suitable  to  the  integration  process  used.  Two  new  variables  and  Vy  are  defined  by 


iX  - V 
dt  “ '^x 


iX  = V 

dt  '^Y 


Substituting  the  above  expressions  in  equations  (2-1)  results  in 


-dr=-«^x 


-HVY-g(Y) 


The  four  first-order,  differential  equations  (2-4)  and  (2-5)  are  the  desired  equations  with  time  as  the 
independent  variable. 


24 


D 162-1 0026-1 


2.2  Transformation  of  Basic  bciuations 

For  reasons  which  will  become  clear  later  in  the  document,  it  is  desirable  to  transform  the  differential 
equations  of  motion  into  two  other  forms.  One  form  uses  downrange  X as  the  independent  variable 
(X-based),  and  the  other  form  uses  altitude  Y as  the  independent  variable  (Y-based). 

The  terms  X-based,  Y-based,  and  time-based  are  adopted  to  ease  the  terminology  in  the  text. 


The  transformation  to  make  X the  independent  variable  will  be  performed  first.  Multiplying  the 

equations  in  (2-5)  by  or 

dX  Vx 


(Note; 


dt_  _1_ 
dX-  Vx 


). 


dX 


= -H 


(2-6) 


dVy  -11  Vy  -g(Y) 
dX  " Vx 


(2-7) 


Similarly,  the  equations  in  (2-4)  can  be  transformed  to 


(2-8) 


(2-9) 


Equations  (2-6)  through  (2-9)  now  represent  a complete  formulation  of  the  ballistics  equations  in 
terms  of  the  variables  Vx  , Vy  , Y and  t , with  X as  the  independent  variable.  Note  that  the 

equations  in  (2-6),  (2-7),  and  (2-8)  are  uncoupled  from  (2-9),  the  equation  involving  time.  Therefore, 
the  first  three  equations  (2-6),  (2-7),  and  (2-8)  can  be  solved  simultaneously  without  solving  for  the 
variable  t by  equation  (2-9).  Figure  5 shows  the  salient  features  in  the  X-Y  plane. 


and 


dX 


'X 


iL 

dX 
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Y 


Figure  5.  Saliant  Features  in  X-Y  Plant 


One  can  form  a set  of  differential  equations  with  Y as  the  independent  variable  in  a similar  fashion. 
Multiplying  the  equations  in  (2-5)  by  ^ or  ^ 

(Note;  -p7  = -77-  ) results  in 
d Y V Y 


(2-10) 


Also,  from  (2-4) 


(2-11) 


(2-12) 


(2-13) 


D 162-1 0026-1 


Hquations  (2-10)  through  (2-13)  represent  a eomplete  formulation  of  the  ballistics  equations  in  terms 
of  the  variables  , Vy  , X and  t , with  Y as  the  independent  variable.  Note,  as  before,  the 

equation  involving  time  (2-13),  is  uncoupled  from  the  other  three  equations. 

Also  note  that  the  equations  with  X as  the  independent  variable,  (2-6)  through  (2-9),  become 
indeterminate  when  is  zero  (e.g.,  when  the  projectile  is  coming  straight  down),  because 

appears  in  the  denominator.  Similarly,  the  equations  with  Y as  the  independent  variable,  (2-10) 
through  (2-13),  become  indeterminate  when  Vy  is  zero  (e.g.,  at  level  release).  Consequently,  the 

two  sets  of  equations  must  be  used  accordingly. 

The  basic  set  of  differential  equations  in  the  integration  procedure  developed  here  is  X-based.  This 
choice  was  based  on  experimental  evidence  which  showed  that  this  set  of  eijuations  gave  the  most 
favorable  results.  A sample  of  this  experimental  evidence  is  given  in  Section  2.4.  The  other  two  sets 
of  equations  are  used  in  special  eases.  When  time  fuzing  is  involved,  it  is  convenient  to  use  the  time- 
based  equations  for  the  part  of  the  trajectory  between  release  and  fusing.  Y-based  equations  work 
best  when  the  trajectory  is  very  steep;  i.e..  the  projectile  is  coming  almost  straight  down. 

2.3  Integration  Procedure 

As  pointed  out  earlier,  the  differential  equations  describing  a projectile’s  trajectory  are  not  analytically 
integrable  if  an  accurate  model  of  function  II  is  used.  Consequently,  a numerical  integration  pro- 
cedure must  be  used.  The  approach  taken  in  the  generation  of  ballistic  tables  is  to  integrate  the 
differential  equations  numerically,  using  time  as  the  independent  variable.  Many  steps  are  taken  to 
insure  the  desired  accuracy.  (A  typical  step  length  is  0.01  second.)  This  metliod  is  accurate  and 
tlexible,  but  requires  too  much  computer  time  to  be  suitable  for  weapon  delivery  application. 

The  procedure  which  proved  to  be  most  successful  in  eliminating  this  shortcoming  is  to  use  the 
fourth-order,  Runge-Kutta  integration  formulas  with  the  X-based,  differential  equations  to 
integrate  the  trajectory  to  the  vicinity  of  the  impact  point.  Integration  of  the  Y-based  equations 
is  then  used  to  locate  the  impact  point  precisely,  since  its  Y-coordinate  is  known.  An  estimate 
of  the  required  computer  time  and  memory  for  using  this  procedure  to  calculate  a weapon  impact 
point  is  given  in  Section  4. 1 . 

l ime-based  and  Y-based  equations  are  used  with  the  X-based  differential  equations  in  special 
cases.  F'or  example,  it  is  convenient  to  use  the  time-based  equations  for  the  powered  segment 
of  unguided  rocket  trajectory  and  for  the  first  segment  of  the  time-fuzed  dispenser  weapon 
trajectory.  It  is  desirable  to  switch  to  the  Y-based  equations  when  a weapon  trajectory  becomes 
very  steep  like  at  the  end  of  the  long  trajectory  of  a drogued  bomb. 

Details  of  the  integration  procedure  and  the  fourth-order,  Runge-Kutta  formulas  are  described  in 
Sectioi  ■’.3. 1 . Application  of  these  formulas  to  the  X-based,  differential  equations  is  described 
in  Section  2.3.2. 

2.3.1  The  Runge-Kutta  Integration  I'ormulas 

I'he  Runge-Kutta  formulas  provide  a step-by-step  method  of  finding  dependent  variable  values  at 
given  intervals  of  the  independent  variable.  That  is,  values  of  the  dependent  variables  at  a 
particular  value  of  the  independent  variable  are  found  in  terms  of  their  values  at  the  previous  step. 
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Knowing  the  differential  equations  of  a given  system  and  the  initial  conditions,  this  process  may  be 
continued  until  the  desired  range  is  covered.  This  integration  method  requires  no  preliminary 
differentiation  of  the  dependent  variables  beyond  those  prescribed.  This  makes  it  particularly 
useful  if  certain  coefficients  in  the  differential  equations  are  empirical  functions  for  which 
analytical  expressions  are  not  known. 

The  standard,  fourth-order,  Runge-Kutta  formulas,  which  yield  very  good  results  even  for  coarse 
integration  intervals,  are  described  for  a system  of  three  ordinary  differential  equations.  Consider 
the  system  of  differential  equations. 


dY  I, 

dx  = 

(X, 

Y, 

Vx,  Vy) 

dX 

(X, 

Y, 

Vx.  Vy) 

(2-14) 

II 

(X, 

Y, 

Vx,  Vy). 

For  example,  if  the  X-based  equations  describing  a projectile’s  trajectory  were  to  be  solved,  then  the 
above  equations  wouUl  be  of  the  following  form. 


D 162-10026-1 


where 

Kl  - h Fj  (X|^  + -j  h.  Y|^  ^■  Ky,  M^,  Vy^  +-jNq) 

K2  = hF,  (X„+-^h,  K,,  Vx^+  T Vy^+|n,) 

K3  = hF,  (X„  + h,  Y„  -I-  Kl.  Vx  + M^.Vy  + N^) 

" ' " " (2-17) 

and  the  values  of  , M,  , M-)  , and  M3  are  coniputetl  by  exactly  the  same  formulas  as  the  corre- 
spondingly subscripted  K’s , except  that  F-i  is  used  instead  of  F,  . Similarly,  Nq  , N,  , N-j  and 
N3  are  computed  using  F3  instead  of  F,  . 

The  fourth-order,  Rungc'-Kutta  formulas  can  be  expanded  in  a similar  fashion  to  a system  of  four 
first-order,  differential  equations. 

When  applying  numerical  integration  formulas,  it  is  necessary  to  decide  the  size  of  the  integration 
steps.  For  this  problem,  where  very  few  integration  steps  are  taken  (i.e.,  on  the  order  of  five),  the 
step  size  depends  largely  on  the  impact  range  of  the  specific  trajectory.  Therefore,  it  is  necessary  to 
have  an  estimate  of  the  impact  range  before  starting  integration. 

One  way  of  obtaining  this  estimate  is  to  use  the  vacuum  trajectory  in  finding  the  impact  range.  This 
works  well  for  low-drag  bombs;  however,  this  estimate  may  be  off  by  as  much  as  a factor  of  two  or 
three  in  the  high-drag  bomb  case.  Fstimating  the  impact  range  will  not  be  a problem  in  practice. 

Since  the  impact  range  will  be  calculated  repetitively,  one  will  always  have  a good  estimate  of  the 
impact  range  once  the  process  is  started. 

2.3.2  Integration  Algorithm 

Originally,  the  integration  was  performed  in  the  following  manner.  The  fourth-order,  Runge-Kutta 
formulas  were  applied  to  the  X-based,  differential  equations,  (2-6),  (2-7),  and  (2-8),  with  a test 
performed  after  each  integration  step  to  determine  whether  the  bombs  were  below  target  altitude. 

When  the  bomb  was  found  to  be  below  this  altitude,  a switch  was  made  to  the  Y-based  equations. 
(2-10),  (2-1 1 ),  and  (2-12).  Then  a single  integration  step  was  taken  to  the  target  since  the  integration 
step  size  in  the  Y-based  system  was  known  exactly.  (See  Figure  6.)  Using  this  procedure,  it  was 
found  that  the  solution  accuracy  for  drogued  weapons  degraded  rapidly  as  release  altitude  went 
above  2,000  feet.  was  so  close  to  zero  at  impact  that  some  of  the  approximations  internal  to  the 

Runge-Kutta  process  had  either  zero  or  negative.  As  appears  in  the  denominator,  it 

renders  the  expression  indeterminate.  The  problem  was  resolved  by  adding  two  logical  tests  to  the 
calculations. 
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Figure  6.  Equation  Rooting  - Integration  in  Y 

Coordinates  are  switched  if  the  vertieal  distanee  the  bomb  fell  in  the  last  downrange  integration  step 
is  greater  than  the  distance  to  the  ground.  If  the  next  integration  step  would  put  the  bomb  far 
below  the  ground,  this  replaces  a long  step  in  X and  a long  step  in  Y with  a short  step  in  Y. 

Calculation  time  is  saved  and  accuracy  is  improved. 

The  other  test  is  on  the  ratio  of  velocity  components.  When  the  vertical  velocity  is  larger  than  the 
horizontal  velocity,  the  coordinates  can  be  switched  to  advantage.  Different  steepness  criteria, 
ranging  from  1 to  5,  have  been  tried.  Now  the  test  is  of  this  form; 

-Vy 

If  — — > 2,  switch  coordinates. 

Since  this  switch  may  take  place  some  distance  from  the  impact  point,  it  is  sometimes  necessary  to 

take  more  than  one  step  in  Y . _ j 

A very  simple  method  is  used  now.  If  the  estimated  range  was  divided  into  steps  and  “i”  is  the 

number  of  steps  taken  so  far,  then  the  number  of  Y steps  to  take  is:  I ^ 

r r 

My  = Nj^  - i + 1.  if  (Nj^  - i ) > 0,  or  . j 

My  ~ if  (^x  ■ i)  ^ 0 • T’ 

; t 

n 

ti 
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The  addition  of  these  two  tests  has  made  it  possible  to  calculate  Snakeye  trajectories  from  releases 
as  high  as  70,000  feet  without  encountering  singularities. 

One  final  topic  concernuig  the  integration  algorithm  needs  to  be  covered.  This  concerns  the  question 
of  how  many  integration  steps  are  required  in  the  process  of  solving  for  the  weapon  impact  point. 
The  final  result  is  more  accurate  with  a large  number  of  integration  steps,  but  computer  time  is 
increased.  The  final  results  may  not  be  sufficiently  accurate  if  too  few  steps  are  taken.  Therefore, 
when  choosing  the  number  of  integration  steps  to  be  used,  a trade-off  exists  between  accuracy  and 
computation  time. 


It  was  decided  on  the  basis  of  experimental  results  to  use  five  integration  steps  now.  A typical 
example  of  these  results  is  given  in  Figure  7.  This  number  of  steps  is  somewhat  arbitary  and  could 
be  adjusted  easily  for  specific  applications. 
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Figure  7.  Integration  Convergence 

2.4  Comparative  Efficiencies  of  Coordinate  Systems 

The  discovery  leading  to  the  development  presented  in  this  document  was  that  the  Taylor’s  series 
with  downrange  as  the  independent  variable  (downrange-base)  converged  faster  than  that  with  time 
as  the  independent  variable  (time-base).  However,  no  work  was  done  to  demonstrate  directly  that 
trajeci  '■y  integration  converged  more  rapidly  with  downrange  as  the  independent  variable.  Since 
there  are  advantages  to  using  time  as  the  independent  variable,  the  comparative  efficiency  of  the 
two  approaches  was  investigated. 
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A computer  program  was  devised  to  calculate  weapon  impact  points  with  the  X-based  and  time-based 
equations  using  identical,  fourth-order,  Runge-Kutta,  integration  formulas.  Pains  were  taken  to 
ensure  exact  comparability  of  results  to  cases  using  the  standard  algorithm.  The  Mk  82S  drogued 
bomb  was  the  primary  weapon  tested.  A comparison  of  results  is  given  in  Figure  8.  The  magnitude 
of  the  error  in  the  computed  range  is  plotted  for  a given  number  of  integration  steps.  The  two  solid 
lines  compare  the  errors  for  release  at  a 500-foot  altitude;  the  dotted  lines  are  fur  release  at  a 3,000- 
foot  altitude.  Note  that  the  errors  are  plotted  on  a logarithmic  scale.  Time-based  integration  for 
this  weapon  requires  almost  three  times  as  many  integration  steps  to  achieve  the  same  accuracy. 

This  comparison  holds  over  a wider  range  of  release  conditions. 

Two  effects  contribute  to  the  difference  in  the  case  of  high  drag  bombs:  ( 1 ) the  bomb’s  deceleration 
is  more  nearly  constant  when  considered  as  a function  of  downrange,  making  it  easier  to  approximate 
its  average  over  an  integration  step;  (2)  due  to  comparative  placement  of  the  integration  steps,  the 
bomb  is  going  much  faster  at  release  than  at  impact  and  travels  much  more  than  one  n-th  of  its  trajec- 
tory in  the  first  one  n-th  of  its  time  of  fall.  The  first  part  of  the  trajectory  also  displays  the  highest  ^ 

deceleration  rates.  Therefore,  the  time-based  integration  takes  longer  steps  in  the  more  critical  portion  j 

of  the  trajectory  than  does  downrange-based  integration  with  the  same  number  of  steps.  This  factor 
in  the  accuracy  difference  could  be  moderated  greatly  with  a suitable  method  of  picking  unequal 
length  time  steps.  However,  a number  of  complexities  arises  with  the  use  of  unequal  steps  and  this  ! 

approach  has  not  been  pursued. 

Also  some  comparative  evaluation  was  done  for  the  Ml  1 7 streamlined  bomb.  This  weapon  slows  i 

down  very  little  during  fall,  so  the  forces  are  more  nearly  constant  and  the  integration  intervals  of  i 

the  two  methods  almost  match.  As  a result,  the  comparative  advantage  of  downrange-based  integra- 
tion is  not  as  great,  although  it  remains  superior  to  time-based  integration.  This  leads  to  a possibility 
that  time-based  integration  may  prove  better  overall  for  some  special  applications,  such  as  stream- 
lined. time-fused  cannisters.  i 

Over  the  entire  spectrum  of  trajectory  calculation  problems,  it  seems  that  downrange-base  integration  i 

has  about  a two  to  one  average  advantage.  Different  methods  have  different  logical  advantages  when  ; 

dealing  with  special  situations  such  as  altitude-fusing,  time-fusing,  loft-delivery,  etc.  To  achieve  the 
optimum  for  specific  weapons  one  may  have  to  use  a combination  of  coordinate  systems. 
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Figure  8.  Convergence  Comparison  Range  - Based  and  Time -Based  Integration 
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3.0  ADAPTATIONS  K)  VARIOUS  WTAPONS  AND  SAMPLE  RESULTS 

This  section  contains  the  mathcmatical/logical  adaptations  necessary  to  calculate  impact  points  for 
various  weapons.  Also,  sample  results  and  t|ualitative  descriptions  of  these  are  given  for  a variety  of 
weapon  types.  The  weapon  types  covered  are;  (I)  streamlined  bombs,  (2)  guns,  (3)  drogued 
bombs,  (4) 'cluster  bombs,  and  (5)  unguided  rockets. 

The  tabular  form  in  which  the  sample  results  are  presented  requires  some  explanation.  There  are 
three  entries  for  each  identified  set  of  release  conditions.  The  first  entry  is  the  “range  error”  in 
feet.  This  is  the  difference  between  the  weapon  range  given  in  the  ballistic  table  and  the  range 
computed  by  the  algorithm.  A negative  sign  on  the  error  means  that  the  computed  range  was 
smaller  than  the  table  range.  The  next  entry  is  “angular  error”  in  milliradians.  This  converts  the 
range  error  into  the  equivalent  angular  error  at  the  release  point.  This  angular  form  is  convenient 
for  comparison  with  other  fire  control  system  errors,  most  of  which  are  basically  angular.  A good 
rule  ol  thumb  is  that  if  the  range  error  is  less  than  lO  feet  or  the  angular  error  is  less  than  I milli- 
radian  the  computation  error  will  not  significantly  degrade  the  overall  system  accuracy.  The  third 
entry  is  the  “time  error”  in  percent.  By  including  the  uncoupled  time  equation  in  the  integration, 
the  time  of  fall  can  be  computed  quite  accurately.  Ability  to  compute  the  time  of  fall  accurately 
will  become  important  when  the  algorithm  is  put  to  an  operational  application.  The  target  for  all 
tables  is  at  sea  level  and  there  is  no  wind;  this  matches  the  conditions  in  the  ballistic  tables. 

3. 1 ♦ Streamlined  Bombs 

Calculation  of  streamlined  bomb  impact  ranges  does  not  require  any  adaptation  of  the  process 
described  in  Section  2.3.2.  Sample  results  comparing  the  impact  ranges  calculated  by  the  algorithm 
to  those  from  Navy  ballistic  tables  are  given  in  fables  VII  through  IX  for  the  Ml  1 7,  Mk  82,  and 
Mk  76  bombs,  respectively.  Table  X gives  results  outside  the  normal  A-7E  operating  region  for  the 
Mk  76  bomb.  I he  tables  show  that  the  algorithm  yields  good  results  over  a large  range  of  release 
conditio’is. 

3.2  duns 

The  approach  taken  with  guns  is  to  treat  them  as  streamlined  bombs,  handling  the  muzzle  velocity 
as  a large  ejection  velocity.  Sample  results  for  the  M6l  gun  are  given  in  Table  XI.  The  calculation 
errors  are  defined  as  they  were  defined  for  bombs  in  the  previous  section.  These  sample  calculations 
* show  that  the  bombing  algorithm  produces  excellent  results  without  any  added  logic  or  calculation. 

3.3  Drogued  Bombs 

One  of  the  major  successes  of  the  algorithm  was  successful  calculation  of  the  impact  point  of  a 
drogued  bomb.  These  bombs  are  dropped  in  a streamlined  configuration  and  then  vanes  or  a para- 
chute is  deployed  to  greatly  increase  their  drag  and  retard  their  forward  motion.  Previous  methods 
were  marginal  in  working  with  high-drag  projectiles  and  were  unable  to  cope  with  the  change  in  the 
drag  function.  Presently,  drogued  bombs  are  handled  by  splitting  the  trajectory  into  two  parts, 
using  a simple  approximation  for  the  streamlined  portion  and  the  algorithm  for  the  drogued  portion. 
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Tableau  Ml  17  Bomb  Computation  Accuracy 


RELEASE 

ALTITUDE 

(FT) 

RELEASE  SPEED  (KNOTS) 

200 

400 

600 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

(%) 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

('*») 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

(%) 

1,000 

-t-1 

0.1 

0.0 

0 

0.0 

0.1 

0 

0.0 

0.0 

3,000 

0 

0.0 

0.0 

-FI 

0.0 

0.0 

0 

0.0 

0.0 

5,000 

-FI 

0.1 

0.0 

-FI 

0.0 

0.0 

-1 

0.0 

0.0 

7,000 

-FI 

.0.1 

0.0 

-F2 

0.1 

0.0 

-2 

0.0 

0.0 

9,000 

-FI 

0.1 

0.0 

-FI 

0.0 

0.0 

-FI 

0.0 

0.0 

11.000 

-FI 

0.1 

0.0 

-F2 

0.1 

0.0 

-2 

0.0 

0.0 

13,000 

-FI 

0.1 

0.0 

+ 1 

0.0 

0.0 

-4 

0.1 

0.0 

15.000 

-FI 

0.0 

0.0 

-F2 

0.0 

0.0 

-F2 

0.0 

0.0 

REFERENCE:  NAVAL  WEAPONS  LABORATORY  BALLISTIC  TABLE  NUMBER  101,  JULY  1967 


CONDITIONS:  LEVEL  RELEASE 
NO  WIND 

TARGET  AT  SEA  LEVEL 
NO  EJECTION  VELOCITY 


Table  VIII  Mk  82  Bomb  Computation  Accuracy 


RELEASE 

ALTITUDE 

(FT) 

RELEASE  SPEED  (KNOTS) 

200 

400 

600 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

(%) 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

(%) 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

(%) 

1,000 

0 

0.0 

0.0 

-FI 

0.0 

0.0 

-FI 

0.0 

0.0 

3,000 

0 

0.0 

0.0 

-FI 

0.0 

0.0 

0 

0.0 

0.0 

5.000 

-FI 

0.1 

0.0 

-FI 

0.0 

0.0 

0 

0.0 

0.1 

7,000 

0 

0.0 

0.0 

-FI 

0.0 

0.0 

-F2 

0.0 

0.0 

9,000 

0 

0.0 

0.0 

0 

0.0 

0.0 

-FI 

0.0 

0.0 

11,000 

-FI 

0.1 

0.0 

0 

0.0 

0.0 

-FI 

0.0 

0.0 

13,000 

+ 1 

0.1 

0.0 

-FI 

0.0 

0.0 

0 

0.0 

0.0 

15,000 

-FI 

0.0 

0.0 

0 

0.0 

0.0 

-1 

0.0 

0.0 

REFERENCE:  NAVAL  WEAPONS  LABORATORY  BALLISTIC  TABLE  NUMBER  169,  MARCH  1969 
CONDITIONS:  LEVEL  RELEASE 

TARGET  AT  SEA  LEVEL 
NO  WIND 

NO  EJECTION  VELOCITY 
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' Table  IX  Mk  76  Bomb  Computation  Accuracy 


RELEASE  SPEED  (KNOTS)  | 

RELEASE 

200 

400 

600  1 

ALTITUDE 

RANGE 

ANGULAR 

TIME 

RANGE 

ANGULAR 

TIME 

RANGE 

ANGULAR 

TIME  1 

(FT) 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

(FT) 

(MR) 

(%) 

fFT) 

(MR) 

(%) 

(FT) 

(MR) 

(%) 

1.000 

0 

0.0 

0.0 

0 

0.0 

0.0 

1 

0.0 

0.0 

3.000 

+1 

0.0 

0.0 

-1 

0.0 

0.0 

-6 

0.1 

0.1 

5.000 

-H 

0.0 

0.0 

-2 

0.1 

0.1 

-11 

0.2 

0.1 

7.000 

0 

0.0 

0.0 

-2 

0.1 

0.0 

-15 

0.3 

0.0 

9,000 

+1 

0.1 

0.0 

-3 

0.1 

0.0 

-1-1 

0.0 

0.0 

11,000 

0 

0.0 

0.0 

-4 

0.1 

0.0 

-26 

0.5 

0.0 

13,000 

+ 1 

0.0 

0.0 

-6 

0.2 

0.0 

-34 

0.7 

0.1 

15,000 

0 

0.0 

0.0 

-9 

0.3 

0.1 

-39 

0.8 

0.1 

REFERENCE;  NAVAL  WEAPONS  LABORATORY  BALLISTIC  TABLE  NUMBER  086.  OCTOBER  1%6 
CONDITIONS:  LEVEL  RELEASE 


NOWIND 

TARGET  AT  SEA  LEVEL 
NO  EJECTION  VELOCITY 


Table  X Mk  76  Bomb  Computation  Accuracy  In  Extended  Delivery  Regime 


RELEASE 

ALTITUDE 

(FT) 

RELEASE  SPEED  (KNOTS) 

600 

900 

1,200 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

("») 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

(%) 

10.000 

-22 

0.5 

0.0 

+ 4 

0.1 

0.1 

-152 

2.2 

0.2 

20,000 

-42 

0.8 

0.0 

-43 

0.7 

0.1 

-91 

1.2 

0.3 

30,000 

-48 

0.8 

0.0 

-82 

1.1 

0.0 

-27 

0.3 

0.0 

40.000 

-15 

0.2 

0.2 

-11 

0.1 

0.3 

-39 

0.4 

0.3 

50,000 

-182 

2.2 

3.1 

-123 

1.3 

1.9 

-73 

0.6 

0.1 

60,000 

-38 

0.4 

0.2 

+ 28 

0.2 

0.2 

-36 

0.3 

1.1 

70,000 

-79 

0.7 

4.8 

-17 

0,1 

5.3 

-21 

0.1 

1.8 

REFERENCE;  NAVAL  WEAPONS  LABORATORY  BALLISTIC  TABLE  NUMBER  086,  OCTOBER  1966 
CONDITIONS:  LEVEL  RELEASE 


NO  WIND 

TARGET  AT  SEA  LEVEL 
NO  UECTION  VELOCITY 
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Tabl«  XI  M6/  Cun  Wifh  M56  Round  Computation  Accuracy 


DIVE 

ANGLE 

(DEG) 

RELEASE 

ALTITUDE 

(FT) 

RELEASE  SPEED  (KNOTS) 

300 

450 

GOO 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

(%) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

(%) 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

(%) 

15 

1,000 

-1 

0.1 

0.G» 

-1 

0.1 

0.0 

-1 

0.1 

0.0 

15 

1,250 

-1 

0.1 

0.5* 

-2 

0.1 

0.5* 

-1 

0.1 

0.5* 

15 

1,500 

-2 

0.1 

0.3* 

-2 

0.1 

0.8 

-2 

0.1 

0.4* 

20 

1,000 

mom 

0.0 

3.0 

0 

0.0 

0.0 

0 

0.0 

1.1* 

20 

1,250 

0 

0.0 

0.0 

0 

0.0 

0.7* 

0 

0.0 

0.0 

20 

1,500 

-1 

0.1 

0.5* 

0 

0.0 

0.6* 

-1 

0.1 

0.7* 

20 

1,750 

-1 

0.1 

0.4* 

-1 

0.1 

0.0 

-1 

0.1 

0.5 

‘TIME  ERROR  OF  0.1  SECOND 

REFERENCE:  ARMAMENT  MEMORANDUM  REPORT  64-5,  FEBRUARY  1%4 
CONDITIONS:  TARGET  AT  SEA  LEVEL 


NO  WIND 

3,300  FT  SEC  MUZZLE  VELOCITY 


I As  noted  before,  the  treatment  of  retarded  weapons  is  complicated  by  the  discontinuity  in  the  drag 

1 coefficient  at  drogue  deployment.  If  it  is  assumed  that  the  deployment  occurs  instantaneously  at  a 

time  tf)  which  is  known  or  can  be  calculated  in  advance,  the  differential  equations  (2-1 ) have 
discontinuous  coefficients  at  t^  due  to  the  change  in  the  ballistic  drag  coefficient  . Corre- 
sponding to  t^3  is  point  Xq  at  which  (2-6)  and  (2-7)  have  discontinuous  coefficients.  In  this 

case,  it  is  best  to  split  the  integration  domain  into  two  parts  corresponding  to  the  different  drag 
regions.  (See  Figure  9.) 

If  the  drogue  deployment  occurs  shortly  after  release  as  for  the  Mk  82  Snakeye  drogued  bomb,  this 
difficulty  can  be  overcome  easily.  Since  Xq  's  very  small  compared  to  X|^  and  the  bomb  has  very 

little  drag  in  the  first  region,  it  is  sufficient  to  integrate  from  release  to  Xq  in  one  step.  Normally, 
tj)  is  given  as  a function  of  the  dynamic  air  pressure.  The  point  Xq  and  the  values  of  , Vyj) , 
and  Yq  must  be  computed  before  the  Runge-Kutta  procedure  can  be  used  to  integrate  from  Xq  to 

Xj^ . It  is  probable  that  due  to  shortness  of  tp  (0.5  second)  the  vacuum  solution  would  suffice.  How- 
ever, it  is  not  difficult  to  account  for  most  of  the  drag  effects  in  the  interval  (0,  Xp)  as  follows. 

r ‘ 
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X p (DROGUE  DEPLOYMENT)  X ^ (IMPACT) 


Figun  9.  Drag  Discontinuity 


Xj)  is  expanded  first,  in  a power  series  in  t about  the  point  t = U . That  is; 


X(t)=X(0)  + X(())t  t-+... 


(3-1) 


The  coefficients  are  computed  readily  from  the  release  conditions.  Denoting  the  release  conditions 
by  the  subscript  zero. 


X(0)  = 0,  X(0)=  Vxo-  X(0)  = -HqX(0)  =-HoVxo  (3-2) 
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Retaining  three  terms  in  the  series  and  evaluating  at  tj^  yields 

A similar  expansion  yields  Yq  , 

^ ^0  *D  ^'^XO  --StDCHoVYo  + gCYo)).  (3-4) 

Differentiating  (3-1)  with  respect  to  t and  evaluating  at  tp  gives 

^XD  = '^XO^*  ■ ^0*D^- 

Similarly,  Vyj)  is  given  by 

'^YD  ^ '^YO  ■ *D  ^^o'^YO  ' 

The  initial  value  of  Xj^  , Yj^  , , and  Vyi)  needed  for  the  integration  process  over  the  high-  I 

drag  region  is  given  by  (3-3)  through  (3-6),  respectively.  The  integration  over  this  region  is  performed  i 

as  described  in  Section  2.0. 

Sample  results  for  the  drogued  Mk  82S  bomb  arc  given  in  Table  XII. 

3.4  Cluster  Bombs 

A cluster  bomb,  which  is  a dispenser  containing  many  individual  bombs,  is  generally  released  in  a 
streamlined  shape.  Then,  at  some  predetermined  condition  (altitude,  time  after  release,  slant  range, 
etc.),  the  bombs  are  released  from  the  dispenser.  The  bombs  have  a much  higher  drag-to-mass  ratio 
than  the  dispenser  and  slow  down  rapidly.  Ballistics  agencies  have  developed  drag  tables  for  a 
fictitious  bomb  which  simulates  the  characteristics  of  the  center  of  the  pattern  as  though  it  were  a 
single  bomb.  r - 

r? 
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Tabl»  XII  Mk  82S  Bomb  Computation  Accuracy 


RELEASE  SPEED  (KNOTS) 

RELEASE 

200 

400 

600 

ML  1 1 1 UUl. 

(FT) 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

Cs) 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

(%) 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

(%) 

1,000 

+ 1 

0.3 

0.0 

-FI 

0.2 

0.0 

0 

0.0 

0.0 

2,000 

-FI 

0.3 

0.0 

0 

0.0 

0.0 

0 

0.0 

0.0 

3,000 

0 

0.0 

0.0 

0 

0.0 

0.1 

-1 

0.1 

0.0 

4,000 

-FI 

0.2 

0.0 

0 

0.0 

0.3 

-1 

0.1 

0.1 

5,000 

0 

0.0 

0.0 

-1 

0.1 

0.0 

-3 

0.3 

o.c 

REFERENCE;  NAVAL  WEAPONS  LABORATORY  BALLISTIC  TABLE  NUMBER  010,  OCTOBER  19S4 
CONDITIONS:  20°DIVE 


TARGET  AT  SEA  LEVEL 

NO  WIND 

NO  EJECTION  VELOCITY 

Like  drogued  bombs,  the  cluster  bombs  are  modeled  by  breaking  the  trajectory  into  two  parts.  The 
first  describes  the  trajectory  of  the  streamlined  dispenser  and  the  second  describes  the  trajectory  of 
the  center  of  the  pattern  of  bombs  released  by  the  dispenser.  The  simple  approximation  used  for 
the  first  trajectory  segment  of  the  drogued  bombs  is  not  adequate  for  cluster  bombs  because  the 
trajectory  segment  is  much  longer;  therefore,  a regular  integration  process  must  be  used.  The  sample 
results  given  in  Table  XIII  for  the  CBU-24B/B  weapon  are  quite  good,  considering  that  the  pattern 
is  about  1 ,000  feet  in  diameter. 

Two  types  of  fuzing  systems  are  considered  here:  one  is  altitude  sensitive,  the  other  is  time 
sensitive.  The  basic  coordinate  system  of  the  algorithm  (i.c.,  the  X-based  system)  works  very  well 
for  calculating  both  segments  of  the  trajectory  of  the  altitude-fuzed  weapon.  During  the  first 
segment  of  the  trajectory  for  this  case,  a hypothetical  target  is  assumed  at  the  fuzing  altitude  and  the 
calculations  arc  performed  conventionally.  For  Table  XIII,  two  range  steps  and  one  altitude  integra- 
tion step  were  used.  Next,  the  same  integration  process  was  repeated  for  the  second  segment  of  the 
trajectory  which  is  between  the  fuzing  point  and  the  actual  target.  For  Table  XIll,  three  range  steps 
and  one  altitude  step  were  used.  In  other  words,  two  conventional  trajectory  calculations  are 
performed  to  find  a single  impact  point  of  an  altitude-fuzed  cluster  bomb. 

The  same  method  cannot  be  used  for  time-fuzed  cluster  bombs.  Normally,  the  time  of  fall  is  not 
calculated  when  the  X-based  equations  arc  used.  The  approach  taken  is  to  use  the  time-based 
differential  equations  for  the  first  segment  of  the  trajectory.  This  is  a simpler  alternative  than  an 
elaborate  rooting  scheme.  The  first  portion  of  the  trajectory  can  be  integrated  accurately  because  it 
is  short  and  low  drag.  The  X-based  equations  are  used  for  the  second  segment  of  the  trajectory. 
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Tabim  XIII  CBU-24  S/S  Clusttr  Bomb  Computation  Accuracy 


RELEASE  SPEED  (KNOTS) 

RELEASE 

ALTITUDE 

(FT) 

420* 

500 

600 

RANGE 

ANGULAR 

WBSM 

RANGE 

ANGULAR 

TIME 

RANGE 

ANGULAR 

TIME 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

ERROR 

(FT1 

(MR) 

■HI 

(FT) 

(MR) 

(%) 

(FT) 

(MR) 

(%) 

5,000 

-2 

0,2 

0.3 

-7 

0.6 

0.5 

-33 

2.8 

1.5 

6,000 

-2 

0,2 

0.3 

-8 

0.6 

0.4 

-28 

2.0 

1.3 

7,000 

-2 

0,1 

0.2 

-9 

0.6 

0.5 

-25 

1.5 

1.1 

■11 

-3 

0.2 

0.3 

-9 

0.5 

0.5 

-22 

1.2 

1.0 

9,000 

-4 

0.2 

0.2 

-11 

0.6 

0.6 

-21 

1.0 

1.0 

10,000 

-5 

0.3 

0.3 

-11 

0.5 

0.6 

-19 

0.9 

0.9 

•400-KNOT  SPEED  NOT  IN  TABLE  FOR  HIGHER  ALTITUDES 

REFERENCE:  F-105D  BALLISTIC  TABLES  CBU-24  B B OR  CBU  29  B 'B,  FEBRUARY  1968 

CONDITIONS:  30°OIVE 


TARGET  AT  SEA  LEVEL 
NO  WIND 

6 FT  SEC  EJECTION  VELOCITY 
2,500  FT  FUSE  ALTITUDE 


3.5  Unnuided  Rockets 

The  unguided  rocket  is  more  difficult  because  of  the  added  problems  of  thrust  and  change  of  mass. 
These  functions  are  provided  as  tables  with  time  as  the  independent  variable.  Linear  interpolation  is 
used  for  values  between  points  in  the  tables.  Since  the  primary  variables  are  given  as  a function  of 
time,  it  is  more  convenient  and  accurate  to  use  the  time-based  differential  equations  for  the  powered 
segment  of  the  trajectory.  The  X-based  equations  arc  used  for  the  rest  of  the  trajectory. 

The  necessary  adaptations  for  the  powered  phase  are  covered  next.  At  any  time  tj^  , the  acceleration 
due  to  propulsion  is 

A -Ll 

■'t  - 

Where 

F-j-  = Rocket  thrust  force  at  time  tj^ 

Mj^  = Rocket  mass  at  time  tj^ 
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To  form  a function  comparable  to  the  aerodynamic  drag  function,  H , we  compute 


V 


The  total  drag  function  is 

H =11^+  Hr, 

where  H^  is  the  drag  function  for  unpowered  weapons  (i.e.,  for  unpowercd  weapons  H = H^). 

A program  has  been  written  to  calculate  impact  ranges  for  unguided  rockets.  However,  ballistic  tables 
which  can  be  used  to  compare  results  are  not  available  at  this  time. 

3.6  Sensitivity  to  Integration  Interval  Size 

Results  given  in  the  preceding  sections  are  based  on  dividing  the  exact  trajectory  length  into  the 
appropriate  number  of  integration  steps.  The  argument  was  that  this  is  valid  because  the  trajectory 
length  changes  little  from  one  calculation  to  the  next;  therefore,  use  of  the  previously  calculated 
range  would  yield  almost  the  same  integration  steps.  This  section  deals  with  the  validity  of  this 
assumption.  In  Section  3.6.1,  the  maximum  deviation  of  the  release  conditions  which  might  be 
achieved  by  maneuvers  during  a single  calculation  interval  is  estimated.  The  corresponding  effect  of 
these  deviations  of  the  impact  range  is  found  in  Section  3.6.2.  This  is  done  for  a low-drag  bomb 
(Ml  1 7)  and  a drogued  bomb  (Mk  82S)  for  short,  medium,  and  long  trajectories.  The  maximum 
percentage  impact  range  deviation  due  to  maneuvers  is  estimated  on  the  basis  of  this  data.  Finally, 
these  deviations  are  included  in  determining  the  integration  interval  to  demonstrate  that  they  do 
not  affect  the  results  adversely. 

3.6.1  Maximum  Deviation  in  Release  Conditions  Due  to  Maneuvers 

The  following  assumptions  were  made  in  estimating  maximum  deviation  possible  in  release  conditions 
during  a single  calculation  interval. 

a.  The  A-7  aircraft  is  the  aircraft  of  interest. 

b.  The  calculation  interval  is  0. 1 second. 

c.  Only  positive-load  factors  are  considered,  since  weapons  are  not  released  under  negative- 
load-factor  conditions. 

The  release  conditions  affecting  the  impact  range  are  altitude,  velocity,  and  dive  angle.  Considering 
velocity  first,  it  can  be  reasoned  that  the  change  in  the  magnitude  of  velocity  in  0.1  second  is 
negligible.  Consider  a hypothetical  example  where  an  accelerating  force  equal  to  the  aircraft  weight 
is  applied.  That  is,  this  force  exceeds  that  needed  to  overcome  drag  and  lift.  This  accelerating  force 
would  change  the  velocity  magnitude  in  0.1  second  by 


AV  = a(O.l), 
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I 

I 

i 

where  a is  the  acceleration,  and  is  equal  to  the  gravitational  acceleration  for  the  chosen  example.  i 

Therefore, 


AV  = 3.2  ft /sec 


If  the  average  bomb  release  velocity  is  assumed  to  be  400  knots,  then  the  above  AV  would  be  less 
than  0.5%  of  this  velocity.  Since  the  available  accelerating  force  is  normally  a small  fraction  of  that 
assumed  in  the  hypothetical  example,  it  can  be  reasoned  that  the  corresponding  velocity  deviation 
in  0. 1 second  is  negligible. 

The  maximum  deviation  in  dive  angle  due  to  maneuvers  is  estimated  next.  The  deviation  in  direction 
an  object  is  traveling  during  a time  interval  At  at  a velocity  V and  pulling  a turn  acceleration 
a is. 


where 


A0  = radians, 
r 


Ad  = AtV  and  r = V^/a,  or, 
AO  = ^ At. 


Therefore,  the  maximum  deviation  in  dive  angle  is  achieved  when  the  ratio  a/V  is  the  greatest. 

Looking  at  the  maneuverability  curves  in  the  flight  manual  of  the  A-7A  aircraft,  the  conditions 
which  correspond  to  the  maximum  A6  are  a 7g  load  factor  at  sea  level  with  a corresponding 
velocity  of  385  knots.  These  limits  are  for  a 25,000-pound  clean  aircraft  with  a TF  30-P-6  engine 
on  a standard  day.  Based  on  the  above  information,  it  is  assumed  that  5g  is  the  maximum 
acceleration  that  can  be  pulled  by  an  armed  aircraft. 

Estimating  the  maximum  dive  angle  deviation  in  0.1  second  using  5g  for  a and  385  knots  for 
V in  the  expression  for  A6  it  is  found  that 


I 


Ae. 


max 


1.4° 


Next,  the  altitude  change  in  one  calculation  time  interval  is  estimated  by  assuming  that  the  dive 
angle  at  the  beginning  of  this  interval  is  maintained  throughout  the  interval.  Inspecting  the  value 
of  ^©niax  ’ considered  a good  assumption.  Therefore,  knowing  the  aircraft  velocity 

V and  dive  angle  © , the  altitude  change  in  0.1  second  is. 


ay  = (0.1  )V  sin© 


ii 
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Summarizing  the  results  of  the  estimated  maximum  deviations  in  the  release  eonditions  in  one 
ealculation  interval; 

a.  The  deviation  in  the  magnitude  of  velocity  is  negligible. 

b.  The  maximum  deviation  in  the  dive  angle  is  1.4  degrees. 

c.  The  ehange  in  altitude  is  AY  = (0.1)  V sin  0. 

3.6.2  hftects  of  Maneuvers  on  Impact  Range  Calculations 

Recall  that  the  objective  here  is  to  demonstrate  the  usability  of  the  impact  range  from  the  previous 
calculation  in  obtaining  the  integration  step  size  for  the  present  calculation.  This  is  done  by 
demonstrating  that  the  results  of  imi^act  range  calculation  obtained  in  this  manner  are,  for  all  prac- 
tical purposes,  the  .same  as  those  obtained  when  the  exact  impact  range  is  used  to  obtain  the 
integration  step  size. 

I he  maximum  deviation  in  release  conditions  due  to  maneuvers  was  estimated  in  the  previous 
subsection.  Next,  the  effects  of  these  deviations  on  the  impact  range  calculations  are  presented  in 
two  steps:  ( 1 ) the  effect  of  these  release  condition  deviations  on  the  impact  range  is  estimated 
to  determine  how  much  the  impact  range  might  deviate  in  one  calculation  interval;  (2)  the  true 
impact  range  is  perturbed  by  this  amount  and  is  used  to  obtain  the  integration  interval  size. 

Maximum  impact  range  deviation  in  one  calculation  interval  is  estimated  for  a low-drag  bomb  (.Ml  17) 
and  a drogued  bomb  (Mk  82S)  for  short,  medium,  and  long  trajectories.  Results  for  these  cases  are 
given  in  fable  XIV.  In  obtaining  these  results,  the  effects  on  the  impact  range  due  to  altitude  and 
dive  angle  deviations  are  assumed  uncoupled,  since  the  coupling  produces  only  second-order  effects. 
As  mentioneil  previously,  only  positive-load  factors  are  considered.  On  the  basis  of  the  results  in 
fable  XIV,  a conservative  5'7  is  assumed  across  the  board  as  the  maximum  possible  deviation  in  the 
impact  range  due  to  maneuvers.  Note,  for  long  trajectories  the  maximum  deviation  is  more  like  I'f . 

Finally,  fables  X\'  and  XVI  show  calculation  errors  in  impact  ranges  for  the  Ml  17  and  Mk  S2S 
bombs  using  integration  steps  baseil  on  perturbed  impact  ranges  erjual  to  the  true  impact  range  minus 
S'r  of  this  range  to  account  for  maneuver  deviations.  1 he  points  in  these  tables  correspond  to  the 
points  in  fab’es  VII  and  VIII.  which  were  obtained  by  using  the  exact  impact  range  to  establish 
the  inlegratio  , intervals,  fhese  results  show  that  even  under  the  most  extreme  maneuvers  the 
errors  in  impact-range  calcidations  caused  by  those  maneuvers  will  be  acceptable. 

3.7  Advanced  Applic;itions 

Beyond  the  application  to  an  aircraft  fire  control  system,  use  of  the  algorithm  can  provide  improve- 
ments in  several  areas,  l our  areas  are  listed  below: 

3.7.1  Ballistic  fable  Calculation 

Current  staiularils  of  accuracy  coulil  be  maintained,  with  ilecreasing  computer  costs,  if  the  algorithm 
ajiproach  to  trajectory  integration  were  substituted  for  the  methods  in  use.  However,  as  ballistic 
tables  are  not  calculated  often,  the  cost  saving  may  not  justify  the  cost  of  developing  new  computer 
progr.ims. 
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Tabit  XIV  Dtyiation  in  Imppct  Range  Due  to  Altitude  Deviation  in  0.1  Second  and  a 1.4°  Deviation 


in  Dive  Angle 

RELEASE  CONDITIONS 

ESTIMATED  MAXIMUM 
RANGE  DEVIATION 

DUE  TO; 

TOTAL 
DEVIATION 
DUE  TO 
MANEUVERS 
(%) 

BOMB 

TRAJECTORY 

LENGTH 

ALTITUDE 

(FEET) 

VELOCITY 

(KNOTS) 

DIVE  ANGLE 

(DEGREES) 

ALTITUDE 

DEVIATION 

("o) 

DIVE  ANGLE 
DEVIATION 

(%) 

LOW 

SHORT 

1,000 

200 

-50 

-1.4 

-1-  4.4 

+ 3.0 

DRAG 

MEDIUM 

8,000 

400 

0 

-0.0 

+ 2.1 

-1-2.1 

IMIIZI 

LONG 

15,000 

600 

+ 45 

0.0 

+ 1.0 

-Fl.O 

HIGH 

SHORT 

1,000 

200 

-5u 

-1.8 

-f-  n 

-(-2.5 

DRAG 

MEDIUM 

2,500 

400 

-30 

-0.5 

+ 2.2 

-1-1.7 

(Mk  82S' 

LONG 

5,000 

600 

0 

0.0 

-(-  0.9 

-(-0.9 

NOTES;  (0  NEGATIVE  DEVIATIONS  INDICATE  THAT  THE  IMPACT  RANGE  BECAME  SMALLER 

(2'  ONLY  THOSE  DIVE  ANGLE  DEVIATIONS  CAUSED  BY  POSITIVE  LOAD  FACTORS  WERE  CONSIDERED 


Table  XV  M//7  Bomb  Computation  Accuracy  with  Initial  Range  Estimate  5%  Short 


RELEASE 

ALTITUDE 

(FT) 

RELEASE  SPEED  (KNOTS' 

200 

400 

600 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR' 

TIME 

ERROR 

("c) 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

O’..) 

RANGE 

ERROR 

(FT) 

ANGULAR 

ERROR 

(MR) 

TIME 

ERROR 

(%) 

1.000 

-FI 

0.1 

0.0 

-FI 

0.1 

0.0 

0 

0.0 

0.0 

3.000 

0 

0.0 

0.0 

FI 

0.0 

0.0 

0 

0.0 

0.0 

5.000 

-FI 

0.1 

0.0 

FI 

0.0 

0.0 

-1 

0.0 

0.0 

7,000 

-FI 

0.1 

0.0 

F2 

0.1 

0.0 

-1 

0.0 

0.0 

9,000 

-FI 

0.1 

0.0 

FI 

0.0 

0.0 

FI 

0.0 

0.0 

11,000 

-FI 

0.1 

0.0 

F2 

0.1 

0.0 

-3 

0.0 

0.0 

13,000 

-FI 

0.1 

0.0 

FI 

0.0 

0.0 

-4 

0.1 

0.0 

15,000 

■FI 

0.0 

0.0 

F2 

0.0 

0.0 

F2 

0.0 

0.0 

REFERENCE;  NAVAL  YlEAPONS  LABORATORY  BALLISTIC  TABLE  NUMBER  101,  JULY  1967 
CONDITIONS;  LEVEL  RELEASE 
NO  WND 

TARGET  AT  SEA  LEVEL 
NO  EJECTION  VELOCITY 
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Tabu  XVI  Mk  82S  Bomb  Computafion  Accuracy  with  Initial  Range  Estimate  5%  Short 


RELEASE  SPEED  (KNOTS! 

RELEASE 

ALTITUDE 

(FT> 

200 

400 

600 

ANGULAR 

TIME 

RANGE 

ANGULAR 

ANGULAR 

Hiin 

ERROR 

ERROR 

ERROR 

ERROR 

■HrltMJ 

ERROR 

(FP 

(MR! 

(FT) 

(MR) 

(FT' 

(MR) 

■■ 

1.000 

+ 1 

0.3 

0.0 

-FI 

0.2 

0.0 

0 

0.0 

0.0 

2,000 

-FI 

0.3 

0.0 

0 

0.0 

0.0 

0 

0.0 

0.0 

3,000 

0 

0.0 

0.0 

0 

0.0 

0.0 

-7 

0.9 

0.6 

4,000 

+ 1 

0.2 

0.0 

0 

0.0 

0.3 

0 

0.0 

0.1 

5,000 

0 

0.0 

0.0 

0 

0.0 

0.0 

-1 

0.1 

0.0 

REFERENCE:  NAVAL  WEAPONS  LABORATORY  BALLISTIC  TABLE  NUMBER  010,  OCTOBER  1964 
CONDITIONS.  20°  DIVE 


TARGET  AT  SEA  LEVEL 
NO  WIND 

NO  EJECTION  VELOCITY 


3.7.2  Fire  Control  System  Sinuilutor 


[ 


Otie  type  of  work  whieh  reijuires  the  caleiikition  of  many  ballistie  trajectories  is  the  simulation  and 
evaluation  of  a fire  eontrol  system.  Many  eurrent  simulation  programs  are  forced  to  rely  heavily  on 
i|uestionable  statistical  assumptions  about  the  kind  of  errors  which  will  be  experienced  and  how 
they  w ill  combine  in  weapi)ti  ilelivery.  The  “foree”  is  applied  by  the  sheer  eost  of  running  numerous 
trajectories  in  a Monte  Carlo  simulation.  By  reducing  costs,  perhaps  by  a factor  of  100  compared  to 
standard  ballistic  methods,  the  algorithm  ean  economically  provide  more  accurate  and  detailed 
simulation  and  analysis  of  proposcnl  fire  control  systems. 


3.7.3  Ballistie  Missile  Ciuidance 

As  the  requirements  for  missile  accuracy  are  tightened,  along  with  a need  for  greater  tlexibility.  the 
need  for  u.sing  an  explicit  guidance  techni(|ue  based  on  trajectory  integration  increases. 


In  some  applications,  the  new  integration  method  will  provide  the  key  to  obtaining  this  improved 


capability  without  as  great  an  increase  in  guidance  computer  power  as  would  be  retiuired  otherwise. 


3.7.4  Command-duideil  Botnb 

I he  command-guided  bomb  is  a coiKvpt  which  has  been  studied  as  a means  of  improving  the  accuracy 
of  conventional  weapon  delivery  without  greatly  increasing  the  cost  per  pound  of  warhead  delivered. 
The  iilea  is  to  put  a minimal  control  system  and  command  receiver  on  an  otherwise  ordinary  bomb. 

All  expensive  guidance  computers  and  sensors  for  target  acquisition  and  bomb  tracking  are  in  the 
aircralt.  where  their  ex|icnse  may  he  amortized  over  many  delivered  weapons.  The  bombs  are 
released  in  much  the  s;ime  way  as  unguided  bombs.  Then  the  guidance  system  in  the  aircraft  com- 
putes commands  automatically  to  make  small  corrections  to  the  bomb  trajectory  as  the  bomb  falls. 

I his  can  yield  substantially  improved  accuracy,  compared  to  unguided  bombs,  at  a eost  increase  that 
is  small  compareil  to  the  approaches  used  in  guided  missiles. 
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4.0  AUXILIARY  TOPICS 

4. 1 Computer  Rcciuiremcnts 

Based  on  general  judgment  and  automated  estimates  (Reference  1),  it  was  believed  that  the  algorithm 
was  competitive  in  terms  of  computer  memory  and  time  requirements.  Larly  in  1968  (before  the 
new  atmosphere  model  was  developed),  a better  estimate  was  sought.  To  this  end,  a basic  machine 
language  program  was  coded.  To  make  the  job  easier,  a simple  hypothetieal  instruction  repertoire 
was  assumed.  Less  than  20  different  instructions  were  used,  but  the  machine  was  assumed  to  have 
an  index  register  and  lloating  point  arithmetic.  It  was  assumed  for  this  preliminary  estimate  that 
double  precision  arithiretic  was  not  required.  The  entire  repetitive  portion  of  the  algorithm  was 
programmed,  including  ai)  subroutines.  Initial  functions  were  relegated  to  a separate  routine  which 
was  neither  time  critical  nor  coded. 

fable  XVll  lists  the  characteristics  of  the  main  program  (MP)  and  its  eight  subroutines.  MP  uses 
these  subroutines  to  take  N fourth-order,  Runge-Kutta  integration  steps  with  downrange  (X) 
as  the  independent  variable  and  one  fourth-order,  Runge-Kutta  rooting  step  with  altitude  (Y)  as 
the  independent  variable,  fhe  general  layout  of  the  program  was  chosen  to  minimize  the  storage 
requirement  without  unreasonably  penalizing  the  execution  time.  The  subroutines  are: 

Integrate  in  X (INTX) 

Runge-Kutta  Lquations  in  X (RKX) 

Runge-Kutta  I quations  in  Y (RKY) 

Common  Runge-Kutta  Liiuations  (RKC) 

Drag  force  Lquations  (DR AC) 

Atmosphere  Model  (AfMO) 

Ballistic  Coefficient  Lookup  (BAD 
Power  Subroutine  (PWR) 


Table  XV//  Algorithm  Code  Characteristics 


ROUTINE 

STORAGE 

TIMES 

NUMBER  OF  INSTRUCTIONS  BY  TYPE 

(WORDl 

EXECUTED 

LOGIC 

TRANSFER 

ADD 

MULTIPLY 

DIVIDE 

MP 

68 

1 

11 

31 

11 

B 

1 

INTX 

66 

N 

12 

31 

mM 

BB 

3 

RKX 

7L 

4N 

2 

9 

mm 

7 

2 

RKY 

17 

4 

2 

7 

n 

4 

2 

RKC 

8 

3N  + 1 

2 

1 

5 

0 

0 

DRAG 

15 

4N  4-  4 

5 

7 

1 

2 

0 

ATMO 

91 

4N  + 4 

6 

23 

4 

6 

4 

BAL 

18 

4N  + 4 

6 

9 

7 

1 

PWR 

44 

12N  -1- 12 

8 

9 

7 

12 

B 

TOTAL  (N=n 

349 

375 

638 

319 

430 

84 

TOTAL  (N=:51 

349 

1.135 

1,994 

987 

1.314 

756 
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MP  calls  IN TX  N times  for  the  forward  integration,  and  directly  calls  DRAG,  RKY,  and  RkS  once 
eaeh  for  the  rooting  step.  INTX  ealls  DRAG,  RKY,  and  RKS.  DRAG  calls  ATMO  and  BAL  and 
uses  PWR  to  calculate  a square  root.  ATMO  in  turns  uses  PWR  to  evaluate  the  exponential  function 
in  the  atmosphere  model.  Table  XVII  shows  the  resulting  number  of  times  each  subroutine  is  used, 
depending  on  N . the  number  of  forward  integration  steps  to  be  taken.  At  the  bottom  of  the  table 
are  the  totals  for  each  class  of  instructions,  for  one  complete  program  execution,  for  N = 1 or 
N = 5 . By  taking  these  totals  and  the  typical  instruction  execution  times  for  a given  computer, 
the  approximate  execution  time  can  be  calculated  for  that  computer  if  its  instruction  repertoire 
does  not  vary  severely  from  that  hypothesized.  For  example,  these  instruction  times,  typical  of  the 
IBM  4-Pi  Model  1:C  airborne  computer,  were  used: 


Logical: 

Data  Fransfer: 

Floating  Point  Add: 

Floating  Point  Multiply: 

Floating  Point  Divide: 

I'sing  these  speeds  and  the  totals  in  Table  ,X 

If  N = 1 
If  N = 5 


3 microseconds 
2 microseconds 
5 microseconds 
10  microseconds 
14  microseconds 

, total  execution  times  were  calculated: 

0.462  milliseconds 
20.052  milliseconds 


I he  estimate  ol  34‘>  word  storage  does  not  vary  with  N . Because  of  the  highly  iterative  nature  of 
the  calculations,  a small  change  in  the  (To^rum  can  make  a large  change  in  the  total  execution  time. 
In  the  program  as  coded,  half  the  total  time  is  used  up  in  the  atmosphere  calculations.  This  is  an 
area  where  many  changes  are  possible,  and  it  has  been  the  subject  of  considerable  research  this 
year.  (See  Section  4.2. 1 

More  recent,  but  less  detailed,  estimates  indicate  a possible  five  to  one  reduction  in  the  computer 
time  required.  I his  is  still  for  the  simple  but  powerful  hypothetical  computer.  The  required 
execution  time  on  the  .A-?!  fire  control  computer  has  not  been  estimated.  Two  opposing  effects 
will  change  the  estimates: 

I he  computer  in  the  .A-7F.  is  slower  and  less  powerful  than  the  hypothetical 
computer.  1 his  will  tend  to  increase  both  execution  time  and  memory 
requirements. 

! or  an  operational  application,  more  advantage  can  be  taken  of  the  iterative 
nature  of  the  calculation  to  reduce  the  average  execution  time.  The  pro- 
gramming for  this  will  increase  the  memory  reiiuiremcnt,  however. 
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4.2  Derivation  of  New  Atmospheric  Model 

1 he  available  expressions  deseribing  atinospherie  density  are  too  complex  to  be  programmed  elTi- 
eiently  on  an  airborne  computer.  There  are  three  expressions,  one  for  each  of  three  altitude 
regions  between  sea  level  and  104,‘)86  feet,  it  is  necessary  to  evaluate  either  an  exponential 
function  or  a function  to  a nonintegral  power  in  each  of  these  regions.  Therefore,  a search  for  a 
more  appropriate  representation  of  the  atmospheric  density  was  conducted. 


Ik'fore  proceeding,  a brief  description  of  the  source  of  the  atmospheric  standard  is  given.  Orig- 
inally, the  1‘)62  U.S.  standard  atmosphere  was  used  in  the  calculation  of  weapon  impact  ranges. 
(See  Keference  7.)  After  a more  desirable  representation  of  the  1962  atmosphere  was  developed, 
it  was  decided  to  incorporate  the  1%5  revision  into  this  model. 

In  the  first  104,986  feel  of  altitude  this  revision  consists  of  a slight  modification  in  the  expressions 
and  the  boundaries  between  w hich  these  expressions  apply.  The  result  of  these  variations  is  .seen 
as  a deviation  in  about  the  fifth  digit  w hen  both  expressions  are  calculated. 

file  remainder  of  this  Section  consists  of  two  parts.  Section  4.2.1  includes  an  outline  and  the 
results  of  the  search  for  a more  efficient  representation  of  the  atmospheric  density.  I'he  1 962 
LI.S.  standard  atmosphere  was  useil  at  this  time.  Section  4.2.2  consists  of  a more  detailed 
derivation  of  the  series  expansion  that  was  finally  adopted.  In  this  case  the  1965  revision  to  the 
H.S.  standaid  almos|ihere  was  included. 


4.2. 1 Search  for  an  Accurate  Power  Series 

The  method  used  to  find  an  appropriate  power  series  rejiresentation  of  the  atmospheric  density  was 
tr)  try  several  expansions  and  pick  the  best  one.  Three  power  series  approximations  were  obtained 
by  making  use  of  the  Chebyshev,  Legendre,  and  Taylor  expansions. 

The  1962  U.S.  standard  atmospheric  model  was  used  at  this  point  without  the  1965  modifications. 
The  atmos|iheric  density  eipiations  from  this  model  are  given  below.  These  expressions  can  be  used 
also  for  comparison  with  the  modified  expressions  given  in  Section  4.2.2.  Note  that  when  two  or 
more  constants  appeareil  in  these  expressions  they  were  multiplied  out.  Also,  the  expression  for 
temperature  has  been  incorporated. 


The  variables  involved  in  these  atmospheric  density  expressions  are: 


where. 


= atmospheric  density  (Ib/ft'^) 
lljj  = geopotential  altitude  (ft) 

A = geometric  altitude  (ft) 

A (20,855,531) 
g ' A + ^0, 8^, 5, 5.^1 


(4-1) 
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First  Region: 


Secoiul  Region: 


0 $ 


lly  < 36,089.239 

22.03627125 

1288.15  -(.0019812)(llg)) 


(1  -(6.87  55856XHr^)(llg)l-‘’-^-‘'^^‘'' 

(4-2) 


36,089.239  ^ llg  < 65,616.798 

= (.0227i887182)e''^-^0^’^‘^“^^'0‘^<H-36,089.239) 


I'liird  Region:  65.616.798  $ 11„  < 100.000 

cr 

p _ 1.190681437  [I  + 1,4068775  X lO'^^ll -656l6.798)r-'^‘^  '^^-'^*‘^‘^ 
''  1216.65  + .0003048  (11„  -65616. 798)| 

(4-4) 

Note  that  the  upper  limit  of  100,000  feet  was  chosen  rather  arbitrarily.  The  actual  value  given  for 
this  limit  was  104.986.877  feet. 

The  procedure  for  finding  a power  series  represenation  using  the  Chehyshev  e.xpansion  will  he  out- 
lined next.  In  principle  this  procedure  is  the  s;ime  as  that  used  for  the  Legendre  expansion.  A 
detailed  description  ot  the  Legendre  expansion  is  included  in  Section  4.2.2;  therefore,  the  Legendre 
expansion  will  not  he  discussed  here.  Since  the  Taylor  series  expansion  is  fairl>  well  known,  it  will 
not  he  discussed  either. 

In  preparation  for  the  Chehyshev  expansion,  the  attiiospheric  density  expressions  in  each  altitude 
region  need  to  he  transformed  so  that  the  range  of  the  independent  variable  x is.  -1  x ^ 1 . 

Then,  the  object  is  to  expaiul  P^^  . the  atmospheric  density,  in  terms  of  the  Chehyshev  polynomials 
of  the  second  kiiul. 

These  polynomials  are  given  below  : 


Co(x) 

= 1 

I'l  (X) 

= 2x 

U-.  (X) 

= 4X--1 

U3  ( X ) 

= 8x-'*-4x 

l'4(x) 

= I6x'^-12x-+  1 

11-5  (X) 

= 32x‘'  - 32x-^  + 6x 

1 tc. 

These  (unctions are  orthogonal  with  respect  to  the  weight  function  \1  - x-  . 
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The  orthoi>onality  conditions  arc: 


/x/T  -x-  U,,(x)  U,^^(x)  d 


-1 


X = ?r/2  , 111  = n 

= 0 , 111  / n 


riic  expansion  in  terms  of  the  C’licbyslicv  polynomials  is  of  the  form 


= AC(1)  L)y(x)  + At\2)  U,(x)  + An3)  lJ2(x) 


(4-6) 


(4-7) 


where  the  AC's  are  constants  found  by  making  use  of  the  orthogonality  conditions.  For  example,  to 
find  AC(  1 ).  multiply  (4-7)  through  by  Uq  \ \ - \2  and  integrate  over  the  region  (-1 , 1 ).  Because  of 
the  orthogonality  conditions  and  the  fact  that  the  AC’s  are  constants,  the  result  is: 


-1 


\/l  -X-  Un(x 


Uq(x)  dx  = AC(l)ff/2 


1 herefore. 


I 

AC(l)=-i  y P^(x)  Vl  - X-  U()(x) 


ilx 


Similarlv. 


A((2)  = 


T f Vl  - X-  U|(x)  d> 


1 

A(  (ii)  =Y  J ^sv^x)  Vl  - X-  U„.|(x)  d) 


(4-8) 


(4-d) 


(4-10) 


ll.iving  found  as  many  AC's  as  desired,  the  expression  for  in  (4-7)  is  known  eomplctely.  Since 

the  C's  in  this  expression  are  all  polynomials  in  x.  one  can  eollect  like  powers  of  x to  get 

as  a power  series  in  x.  After  transforming  these  expressions  so  that  the  altitude  A becomes 

the  iiulependcnt  variable,  one  obtains  the  desireil  power  scries  expressions  for  P^,  in  terms  of  the 
altitude  A. 


53 


!)I()2-I0020-1 


I ho  rosultini;  coctTicionts  of  the  dift'oront  powers  of  the  geometric  altitude.  A,  obtained  from  the 
different  expansions,  are  given  in  l ahle  XVI II.  In  each  case,  six  terms  were  used  in  the  expansion. 
Table  Xl.X  compares  the  performance  for  tlie  three  different  expansions.  Note  how  mucli  better 
the  series  resulting  from  the  Legendre  expansion  is  than  the  otlier  two  series. 


Table  XVIll  Coefficients  Of  The  Different  Powers  Of  A In  The  Expansion  Of  P 


CO  EF  FI  Cl  ENTS 


C„  FOR  THE  EXPRESSION 
N 


P = 1 V" 


COEFFICIENTS  FROM  THE  COEFFICI ENTS  FROM  THE  COEFFICIENTS  FROM  THE 
CHEBYSHEV  EXPANSION  LEGENDRE  EXPANSION  TAYLOR  EXPANSION 


7.687393941  X 10"" 
2.407296101  X 10~® 
4.815919194  X 10"'^ 
1.474875550  X 10"^^ 
3.544103811  X 10"® 
-3.3839  20382  X 10  "® 


2.358894826  X 10"' 
-1.649  391491  X 10  "^ 
5.375930550  X IC""* 
-9.448623261  X IQ-'^ 
8.578708979  X 10 
-3.1536  7 2168  X 10 


2.660965125  X 10"' 
-1.360  7 2 27  7 8 X 10-5 
2.928061740  X 10*® 
-3.244227337  X 10-'5 
1.827244849  X IC-"" 
-4.1538  5 7 7 7 8 X 10 


7.647500017  X 10"^ 
-2.237784147  X 10 
2.515493058  X 10"" 
-1.319129466  X 10"'® 
2.9865  3 3408  X 10  "^2 
- 1.606578496  X 10  "28 


1.242282588  X 10 
-5.593402375  X 10  "8 
1.156099115  X 10  -'0 
-1.349502937  X 10 
8.702175740  X 10  “21 
-2.4239  2 3389  X 10 


1.169028307  X lO'l 
-4.705422820  X 10-8 
8.125386  209  X lO'H 
-7.409812969  X 10-18 
3.5  23  4 24S24  X 10-21 
-6.919278720  X 10-27 


7.647500000  X 10-2 
-2.237783949  X 10"8 
5.030975273  X 10"" 
-9.157272234  X 10"'® 
5.03486496  7 X 10  "28 
-5.983560285  X 10"28 


*2.271887183  X 10"^ 
-1.088171034  X 10"^ 
5.222455116  X 10~" 
-2.511405471  X 10"" 
1.210099  127  X 10"'^ 
-5.8  4 23  08  6 2 0 X 10  "2'’ 


IN  THIS  REGION  THE  TAYLOR  EXPANSION  IS  ABOUT  A = 36.151.79  7 INSTEAD  OF  A=  OASIN  OTHER  CASES 

Table  XIX  Comparing  Values  of  Calculated  by  the  Different  Expansions  to 
These  Calculated  by  the  Adopted  Standard 
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I 

I 

s 


I 


I 


1 


4.2.2  Derivation  of  Power  Series  from  Legendre  txpansion 

I'he  series  derived  from  the  Legendre  expansion  was  adopted  since  it  produced  by  far  the  most 
accurate  representation  of  the  atmospheric  density.  A detailed  description  of  the  derivation  of  this 
series  was  not  given  in  the  previous  section  and  is  given  here  with  the  DJ65  modifications  to  the 
standard  atmosphere.  Section  4.2.2. 1 includes  the  l‘)62  U.S.  standard  atmospheric  density  with 
l‘)6.S  modifications.  Transformations  necessary  to  make  air  density  expressions  compatible  with  the 
Legendre  polynomial  are  given  in  Section  4. 2. 2.2.  The  Legendre  expansion  is  described  in  Sec- 
tion 4. 2. 2.3.  and  the  results  are  given  in  Section  4. 2. 2.4. 

4.2.2. 1 Pxpressions  for  Atmospheric  Density.  I'he  expressions  for  the  atmospheric  density 
(in  Ib/ft'^).  including  the  I6fa5  modifications,  are  given  below.  First,  let. 


Mj,  = geopotential  altitude  (ft) 
A = geometric  altitude  (ft) 


(4-11) 


where  the  relationship  between  A and  II  is. 


,,  ^ A(  20,85S,.S31) 

g A + 2(),8.S5,.S3I 


(4-12) 


First  Region,  0 11^^  < 36,089 

p 22.03(>44887  1 1 -((>.87.S5856X1(}'<’)(II  )| 

vv  ~ (:88.l5-(.()019812)(llg)j  ^ 

Second  Region.  36,089  ^ Mg  < 6.3,617 

= (0.227190667)  t-’  -4.80634274X10'‘'(Mg -36,089) 

Third  Region.  65,617^  Mg  ^ 104,987 

n _ (I.19()(,943.S6)  |1  +(l.40(i87745XI0-<’)(Mc-6.3.6l7)r-''^  '<’-^'‘^-‘' 

" 1 2 1 6.6.3  ■k0.()00304H)(  I Ig  -6.3.6 18^1 


(4-13) 


(4-14) 


(4-15) 
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4.2. 2. 2 Trimsrorniation  of  Variables.  The  input  to  the  weapon  control  system,  in  practice,  will  be 
the  altitude  A,  not  Hg  . I herefore,  the  independent  variable  Hg  in  the  expressions  for  the 

atmospheric  density  given  above  will  be  transformed  to  A by  the  use  of  (4-1 2).  Also,  in  each 
region  the  expansion  will  be  about  the  lower  boundary  limit  of  the  region. 

I'he  Legendre  polynomials  used  are  orthogonal  in  the  region  -1  ^ x ^ 1 , where  x is  the 
independent  variable  for  these  functions.  Therefore,  it  is  necessary  to  make  a linear  transformation 
of  variables  in  the  expressions  for  P^,  in  each  region  so  becomes  a function  of  x.  where. 

-I  < X $ 1.  After  the  expansion  in  terms  of  the  Legendre  polynomials  is  completed,  one  can 
transform  back  to  the  variable  A.  Lhe  appropriate  transformation  for  each  region  is 

^ ^ ~1  ^''^minj  ■'  (4-16) 

L max  'mmj  L J 

where  and  A,„j„  denote  the  upper  and  1 wer  boundary  limits,  respectively,  of  A in  each 

region. 

4. 2. 2.3  Legendre  l-.xpansion.  Different  numbers  of  terms  were  tried  for  representing  the  atmo- 
sphere in  the  first  region.  It  was  decided  to  use  five  terms  or  a fourth-order  representation.  In  the 
next  two  regions  six  terms  or  fifth-order  representations  were  used. 

Next,  the  object  is  to  expand  P^,(x),-1  s x < 1 , in  terms  of  the  Legendre  polynomials  which 
are  given  below. 
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The  orthogonality  conditions  are: 


(x)  (x)  dx  = 0 


2m  + 1 


, ni  n 


, m = n 


(4-18) 


riie  expansion  in  terms  of  the  Legendre  runctions  is  of  the  form: 


P^,  = IXd)  Pq  (X)  + 1)C'(2)  P,  (X)  + IX'(3)  IS  (x)+...  (4-16) 


where  the  DC's  are  constants  that  are  found  by  making  use  of  the  orthogonality  conditions.  I'or 
example,  to  find  DC(  1 ),  multiple  (4-16)  through  by  Pq(x)  and  integrate  over  the  region  (-1,1).  The 
result  is: 


I 

J P^,  tx)  Pq  (x)  dx  = 2DC(1) 


(4-20) 


riierefore. 


Similarly. 


IX(l)=-lr  j (x)  Pq  (X)  dx 


IX'(2)  = 4 


lX(n+l) 


(x)  P|  (x)  dx 


I - J P^(x)P„(x)dx 


(4-21) 


(4-22) 


DI6M0026-I 


As  mentioned  before,  five  terms  are  used  in  tiie  expansion  for  the  first  region,  and  six  in  the  other 
two  regions.  I'he  eoeffieient.s,  IK'(I)  through  D('{5)  or  1X'(6) , are  found  by  using  (4-22).  Next, 
the  object  is  to  expand  in  powers  of  ( . fhis  is  done  in  two  steps.  First,  the  result  of 

substituting  (4-17)  in  (4-14)  and  eollecting  like  powers  of  X is; 

P^  = IIK'(l)  - l/2  1)C'(3)+|lX'(5))  + X [1X(2)  - 4 W)  + I 


+ X-  l4  1X'(3)  - ^IX'(.S))  -x-^  [ 4 •’<-’(4)  - IK(6)) 

+ x-^  1^1)C(5)  1 + x-*^  l4rW'((>)l  , 


(4-23) 


or. 


P^^,  = FXd)  + I C'(2)  X +FX(3)  X-  + F;C(4)  x-^  -h  F((5)  x*^  -h  FC'(6)  x^  , (4-24) 


respectively.  1 he  last  term  in  the  above  expression  does  not  exist  in  the  expression  for  the 
lirst  region. 

Next,  it  is  desirable  to  get  P^^  as  a lunction  of  (A  - A,^^jp) . where  A^^jp  is  the  lower  limit  in  each 
regioti.  1 he  relationshiji  between  x and  (A  - A,„j„)  is 


X .S  ( A - A|^.|j,^) -I 


(4-25) 


w here  .S  = 2/(  - A|,,j,,) . (.See  (4-1  b).)  The  result  of  substituting  (4-25)  into  (4-24)  and  collect- 
ing like  I'owers  of  (A  - ) is, 

P^^,  = |l:(’(l)  - FX(2)  + l■X'(3)  - F((4)  + FX'(5)  - FX’(6)| 


+ 

' '*^min^  |1X(I)- 

21X(3) 

+ 3IX(4) 

+ 

(A-A,„in)-llX(3)- 

3FX'(4) 

+ (>FX’(5) 

+ 

(A-A,„i,)-MF((4(- 

4FC(5) 

+ U)FX'({i)| 

+ 

(A-Ami„)‘^|FX(5)- 

51X((>)| 

8*^ 

+ 

(A- A,„j„)-‘^  |FX(b)| 

(4-2(>) 


5« 


1 
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or 


LC(1)  + LC(2)  (A-A„^in)  + LC(3)  (A  - + LC(4)  ( A - A,,,;,,)-^ 


L('(5)  (A-A„^i„)'^+  LC(6)  (A-A„,jn)5 


min' 
(4-27) 


respectively.  Note  that  (4-26)  and  (4-27)  represent  the  six-term  expansion.  I'o  obtain  the  five-term 
expan.sion  simply  set  hC(6)  and,  con.sequently,  L('(6)  to  zero. 

fhe  LC"s  are  the  desired  coelTieients.  A computer  program  has  been  written  to  calculate  these 
L( 's. 

4.2. 2. 4 Results.  The  coefficients  of  (A  - A„,ji^)  for  the  calculation  of  the  atmospheric  density 

and  atmospheric  density  ratio  are  give,  .n  Table  XX.  fable  XXI  describes  the  accuracy  of  the  power 
series  representation  of  the  atmospheric  density  expressions,  fhis  accuracy  is  described  by  giving 
the  smallest  number  of  digits  in  agreement  between  the  developed  power  series  and  the  expressions 
approximated. 

When  using  the  CDC  6600  computer,  this  series  can  be  evaluated  approximately  23  times  faster  than 
the  original  expressions.  As  seen  in  fable  XXI.  the  agreement  between  this  series  and  the  expressions 
being  approximated  is  at  least  seven  digits  for  altitudes  below  37.000  feet  and  five  digits  for  alti- 
tudes below  100,000  feet. 


Table  XX  Coefficients  of  (A  - A min)  in  the  Expressions  for  Atmospheric  Density  and  Density  Ratio 


ALTITUDE 

COEFFICIENTS  OF  (A  - A mini  N 

COEFFICIENT  OF  (A  - Amin)  N 

REGION 

N 

FOR  ATMOSPHERIC  DENSITY 

FOR  ATMOSPHERIC  DENSITY 

(FTl 

(LB  CU  FTl 

RATIO 

0 

7.6475577  X 10"^ 

9.9999949  X 10"' 

0 TO 

1 

-2.2377695  X 10"® 

-2.9261216  X 10", 

2 

2.5148808  X 10"“ 

3.28  84740  X 10"'° 

36.15  1.558 

3 

-1.3144741  X 10"^® 

-1.7188146  X 10"', 

4 

2.84  1 3555  X 10"^^ 

3.7153745  X 10"“ 

0 

2.2718915  X 10"2 

2.9707396  X 10"' 

36,151.558 

1 

-1.087  %09X  10"® 

-1.4226246  X 10"® 

TO 

2 

2.6036  5 2 3 X 10"“ 

3.4045523  X 10"'® 

3 

-4.0790415  X 10"'® 

-5.33  37  8 05  X 10"'® 

65.824.100 

4 

4.32073  7 3 X 10"“ 

5.6498234  X 10"^ 

5 

-2.4239199  X 10“^ 

-3. 1695330  X 10'® 

0 

5.4956781  X lO'^ 

7.1861834  X 10"2 

65,824,100 

1 

-2.6989696  X 10"^ 

-3.5291898  X 10"® 

2 

6.7688823  X 10"'^ 

8.8510333  X 10"'' 

TO 

3 

-1. 10995  12  X 10"'® 

-1.4513792  X 10"'® 

105.518.180 

4 

1.176  1356  X 10"“ 

1.5379224  X lO"^® 

5 

-6.0891352  X 10"“ 

-7.%21917  X 10"" 

NOTE:  THE  SYMBOL  A DENOTES  GEOMETRIC  ALTITUDE  AND  A DENOTES  THE  LOWER  ALTITUDE  LIMIT 
FOR  EACH  REGION. 
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Table  XXI  Comparison  of  the  Power  Series  to  the  Expressions  They  Represent 


ALTITUDE 

REGION  (FT', 

COMPARISON  • 

ORDER  OF 

SERIES 

OTO 

36,151.6 

7 

4 

SECOND 

REGION 

36,151.6 

TO  65,824.1 

5 

5 

THIRD 

65,824.1 

TO  100,000 

5 

5 

REGION 

100,000 

TO  105,518.2 

4 

•NOTE  THIS  COLUMN  GIVES  THE  SMALLEST  NUMBER  OF  DIGITS  IN  AGREEMENT  TO  VdTHIN  ONE 
HALF  IN  THE  LAST  DIGIT 
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5.0  CONCLUSIONS 

I ho  technique  deserihed  in  this  doeunient  for  calculating  weapon  impact  points  applies  to  any 
unpowered,  unguided  projectile.  It  also  applies  to  some  powered  and/or  guided  projectiles  which 
have  a substantial  ballistic  trajectory  segment.  Fhe  reason  for  this  wide  applicability  is  that  an 
accurate  representation  of  the  differential  equations  of  motion  is  intcgratcil  directly.  The  efficiency 
of  the  new  algorithm  combined  with  the  high  speed  of  the  new  computers  permits  both  accuracy 
and  tlexibility  in  airborne  fire  control  calculation. 

It  has  been  established  that  the  algorithm  can  predict  the  trajectories  of  many  types  of  projectiles, 
riie  primary  application,  of  cenirse,  is  in  airborne  fire  control  systems.  It  can  be  used  in  weapon 
release  calculations  for  most  air-to-ground  weapons  and  some  air-to-air  weapons.  The  algorithm 
could  also  be  used  in  direction  of  ground-to-ground  or  ground-to-air  guns,  unguided  rockets,  or 
guidance  calculations  for  guided  weapons  which  operate  in  a semi-ballistic  mode. 

I here  are  .souiul  technical  reasons  to  believe  that  results  equal  to  those  shown  in  this  document  can 
be  proiluced  under  similar  conditions  for  any  weapon  with  similar  characteristics. 
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